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FOREWORD 


This research was supported under contract NAS1-10277 with the National 
Aeronautics and Space Administration's Langley Research Center. The technical 

monitor on the contract was H. Harris Hamilton II. The results are published 
in two parts: 

Part I. - Description of Basic Method (NASA CR-111921) 

Part II. - Description of Computer Program (NASA CR-111922) 

Part I contains a detailed description of the theoretical approach and basic 
equations that were used to obtain a solution to the problem. The present 
part (Part II) contains a detailed description of the computer program and 
is intended to serve primarily as a user’s manual. 
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CALCULATION OF INVISCID SURFACE 
STREAMLINES AND HEAT TRANSFER 
ON SHUTTLE TYPE CONFIGURATIONS 


Part II. - Description of Computer Program 


by Fred R. DeJamette 
and 

Michael H. Jones 
North Carolina State University 

INTRODUCTION 


Part I of this report (NASA CR-111921) presents a description of the basic 
method used to calculate laminar, transitional, and turbulent heating rates on 
arbitrary blunt-nosed three-dimensional bodies at angles of attack in hypersonic 
flows. The present report (Part II) gives a description of the computer program 
used in Part I. The main program and each subprogram are described by defining 
the pertinent symbols involved, a detailed flow diagram, and a complete computer 
program listing. 

Input and output parameters are discussed in detail. Listings are given 
for the inputs used in Part I to compute heating rates on (1) a blunted 15° 
half-angle cone at 20° incidence and Mach 10.6, (2) a blunted 70® slab delta 
wing at 10° incidence and Mach 8, and (3) the HL-10 lifting body at 20° 
incidence and Mach 10. In addition, the computer program output for two 
streamlines on the blunted 15° half-angle cone is listed. 


PROGRAM DESCRIPTION 


Descriptions of the main program and each subprogram are given below along 
with the pertinent symbols used in each program* Program inputs and outputs 
are described in subsequent sections. The complete computer program is listed 
in Appendix A, and a detailed flow diagram is given in Appendix B. 


Main Program 


The main program reads part of the input data, calculates the initial data 
for each streamline on the e-circle around the stagnation point, and then cal- 
culates and prints the heating rates and other pertinent data along each stream- 
line independently of the other streamlines * A simple flow diagram for the 
entire computer program is given in figure 1, whereas a detailed flow diagram 
is given in Appendix B. The program variables are listed below, and the corre- 
sponding symbols used in Part I of this report are given if applicable. 


ALP 

ALPD 

B 

BBAR 

BBARS 

BETA 

BETAD 

CAL 

CPSI 

CST 

CTHM 

DDLB 

DDP 

DDX 


angle of attack, ot (radians) 
angle of attack, a (degrees) 

ratio of principal velocity gradients at stagnation point 
pressure. gradient parameter, & 

pressure gradient parameter at stagnation point, 

coordinate normal to streamline on body surface, £ (radians) 

BETA in degrees 

cosine of ALP 

cosine of PSI 

constant to convert units of body coordinates to feet 
constant used to calculate laminar momentum thickness 


h 

_Li 

3<j> 


36 


it 


3x 


DEPS radius of e-circle around stagnation point 
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FIGURE 1, SIMPLIFIED FLOW DIAGRAM 


/ READ 
INPUT 
DATA 

■ ZZZEZ 

COMPUTE BODY GEOM, COEFF. FOR SPLINE, LOCATE 
STAGNATION POINT, COMPUTE STAGNATION POINT PROPERTIES 



/ READ g FOR 
/STREAMLINE TOi 
BE COMPUTED / 


7 

C0Mp UTE INITIAL VALUES ON e-CIRCLE 


COMPUTE BODY ANGLES AND THEIR DERIVATIVES 

~ T 

COMPUTE PROPERTIES AT EDGE OF 
BOUNDARY LAYER, CHECK FOR TRANSITION, 

& COMPUTE HEAT-TRANSFER RATE 



3ZZ1 


INCREMENT DISTANCE 


ALONG STREAMLINE 





DGB 


1 311 
h 36 


DGP 

DGS 

DGX 

DLPH 

DPB 


311 

3<|> 

DT 

DS 

1L 

3x 



i l? (p/p s ) 


DPHB 

DPS 

DPSIS 

DS 

DSB 

DSS 

DTHB 

DUEDPS 

DUEDS 

DXB 

D2PB 

EM 

EMIN 


1 3i 
h 36 

D(p/P s ) 

DS 

M 

DS 

integration step size along a streamline 

1 12 . 

" h 36 

_ £ 2 . 

“ DS 

1 3_8 
h 36 

DU 
e 

Dty 

DU 
e 

DS 

1 3x 
h 36 

13 Fl 3 > , s] 
h 3B [h 36 (p/l V 

Mach number at edge of boundary layer, M^ 
freest ream Mach number, M^ 
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(EM) 2 


same as EMS 


body radius, f 

body radius at stagnation point, f 


p V n 


u cst 

e 


D_>Mk 
DS '"DS ) 

D(HI) 

DS 


y CST 
. e 


for KP - 0, 


for KP ?! 0 


for KP - 0, 


for KP ?! 0 


for KP = 0, 


for KP ?! 0 


for KP ?! 0 


D(HI) 

DS 


for KP ?! 0 


body angle r 

compressible form factor, H 

enthalpy at edge of boundary layer, h (ft 2 /sec 2 ) 
incompressible form factor, 

maximum integration step size 

2 2 

stagnation enthalpy, H^(ft /sec ) 
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HSP 

HTR 

HW 

IMAX 

ISTG 

K 

KBM 

KG 

KMAX 

KP 

KP0 

KS 

KTR 

L 

LI 

L2 


MR 

NAM 

NCS 

NMAX 

NR 

NT 

PE 

PHID 

PIN 

P0P 

P0RV 


QWS (6006) / (HS-HW) , for a perfect gas HSP is the stagnation point 

2 

heat transfer coefficient (BTU/ft -sec-*R) 

transformed form factor, H 

2 9 

wall enthalpy, h^ (ft /sec z ) 
see Description of Input 

indicator variable, see Diagnostic Messages 
see function RUNGE subprogram, K«RUNGE 


\see Description of Input 


counter for integration steps along a streamline 
indicator variable, LI «* 1 for laminar region, Ll * 2 for 
transitional and turbulent regions 

indicator variable, L2 *= 1 for laminar region, L2 B 2 for last 

integration step before the beginning of transition, L2 = 3 for 

interior of transition region, and L2 *= 4 for last integration 

step of transition region 

see function RUNGE subprogram 

see Description of Input 

see subroutine SPCH 

see Description of Input 

number of differential equations to be integrated in RUNGE subprogram 

indicator variable, NT « 1 for laminar region, NT - 2 fully turbulent 

region, NT « 3 for laminar region over entire body 

2 

surface pressure, p(lb/ft ) 

circumferential angle, tp (degrees) 

2 

frees tream pressure, p^ (ib/ft ) 
pressure ratio, p/P s 

P /(p /) 

S 3 00 
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PNEWT 


PS 

PSI 

Q0QS 

Q0QSB 

QWB 

QWL 

QWLB 

QWS 

QWSB 

QWT 

QWTB 

RH0E 

RH0IN 

RH0S 

RH0SW 

RNF 

RNM 

RN1 

S 

SAL 
STR 
v XE 
THE 
THM 
TIN 
TRE 
TRI 
TW 
HE 
U0U 
' VIN 



stagnation pressure, p (lb/ft 2 ) 
angle 

heating rate ratio, q^/q^ 

Q0QS (1 - ZW)/(„85 + olSTE - ZW) , for a perfect gas Q0QSB is the 
ratio of the local to the stagnation point heat transfer coefficient 

(Q0QSB) (QWSB) , (BTU-sac/f t 4 ) 

laminar heating rate, q y (BTU/f t 2 --sec) 

laminar value of QWB 

stagnation point heat-transfer rate, q^ (BTU/f t 2 -sec) 

QWS/ (HS - HW) >S 

turbulent heating rate, q (BTU/f t 2 -sec) 

turbulent value of QWB 

density at edge of boundary layer, p (slug/ft 2 ) 

freestream density, p M (slug/ft 3 ) 

stagnation density, (slug/ft 3 ) 

wall density at stagnation point, p w (slug/ft 3 ) 

Reynolds number per foot, p g U Jv (1/ft) 

mom an turn thickness Reynolds number, Rn 
. m 

integrated unit Reynolds number, Rn^ 

distance along a streamline from stagnation point 
sine (ALP) 

STR*Y(1) for KTR— 1, STR=RNM for KTR-0, or STR=RN1 for KTR=1 
enthalpy ratio, t 

e 

inclination angle of inviscid surface streamline, 0 
laminar momentum thickness, G 

m,lam 

freestream temperature, T 

J- see Description of Input 
wall temperature, *R 

velocity at edge of boundary layer, U (ft/sec) 
velocity ratio, U /V 
freestream velocity, V (ft/sec) 

^ CO f 
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Vise coefficient of viscosity at edge of boundary layer, (slug/f t-sec) 

VISS coefficient of viscosity at stagnation point, y g (slug/f t-sec) 

VISW wall coefficient of viscosity, (slug/f t-sec) 

WFN weighting function, w f 

X0 axial location of stagnation point, x g 

YY y-coordinate 

Y (I) q / S F(I)DS, where I = 1,2,*.., NR 

2W wall enthalpy ratio, S 

ZZ z-coordinate 


Function RUNGE Subprogram 

Function RUNGE is a subprogram which uses the fourth-order Runge-Kutta 
method to integrate a system of NR first-order ordinary differential equations 
with a variable step size. As a criterion for varying the computing interval, 
the differential equations are integrated over an interval of step si2e H first 
and then over the same interval with two step sizes of H/2. The two solutions 
are then compared to give an estimate of the error for each variable. If any 
error is larger than EPS = 10 \ these answers are discarded and the computing 
interval H is halvedn If all of the error estimates are less than EPS, the 
answers are allowed and the integration process continues. In addition, the 
step size is either doubled or set equal to HMAX, whichever is the smaller, 
for the next integration cycle. 

The function RUNGE is used in the Main Program and has the form 
K-RUNGE (Y , F , S , DS ,NR ,HMAX ,MR) 

where . K. can have one of three values : 

K « 0 implies the completion of an integration cycle 

K « 1 implies the integration cycle has not been completed 

—8 

K » 2 implies the step size has been halved down to a value below 10 
The dummy arguments used in the subprogram itself are 

FUNCTI0N RUNGE(Y,F,X,H,N,HMAX,MR) 
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where 

Y array of dependent variables to be integrated 

F array of first derivatives of the dependent variables, F “ — 

X independent variable (distance along a streamline) 

H integration step size 

N number or differential equations 

HMAX maximum integration step size 

MR indicator variable, if MR = 1 the previous integration interval is 
recomputed with e new step size H determined in the main program 

Subroutine STAGN 

Subroutine STAGN locates the Newtonian stagnation point and calculates the 
two principal radii of curvature there. A call to STAGN has the form 

CALL STAGN(ALP > IMAX,KMAX,X0,RPER,B,ISTG) 

where the arguments are 

ALP angle of attack, a (radians) 

1MAX “I 

ksee Description of Input 
axial location of stagnation point 

RPER principal radius of curvature in a plane perpendicular to plane 

of symmetry at stagnation point, R^, 

B ratio of principal radii of curvature, R^/Rn 

ISTG indicator variable, -1 for negative angle of attack, 0 for stagnation 

point not found, or +1 for stagnation point found 

Other program variables are 
FF body radius, f 

RlAR principal radius of curvature in plane of symmetry at stagnation 

point, R n 


Subroutine NEWP 

This subroutine calculates the modified Newtonian pressure and its first 
and second derivatives with respect to both the axial (x) and circumferential 
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(<jO directions* A call to NEWP has the form 


CALL NEWP(AA,F) 


where the dummy arguments are 


AA same as PNEWT in main program, AA 83 1 - 

F body radius, f 

Other program variables are 
CA cosine (ALP) 

CP cos <j> 

2 

FF body radius squared, f 


3x 2 34> 


ratio of static to stagnation pressure, p/p c 
3(p/P_) 

o 


3 <P/P S ) 


3(p/P s ) 


3 (P/P s ) 


3 (p/Pg) 
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sine (ALP) 
sin (j) 


Variables in labeled C0MM0N /C0M 2/ and /C0M 5/ are calculated in subrou- 
tine STRMLN. 


Subroutine STRMLN 

For a given axial (x«XX) and circumferential 4=PPH) location on the 
body, this subroutine calculates the body radius (f ) , body angles T and 6 
and their derivatives in the axial and circumferential directions. A call 
to subroutine STRMLN has the form 

CALL STRMLN (XX , PPH , DLPH , DDP , DDX , GM , DGP , DGX , F , IMAX , KMAX) 


where 

■ the arguments are 

XX 

axial position, x 

PPH 

circumferential position, $ 

DLPH 

body angle 5 

DDP 

35, * 

— £ 

3$ 

36 y 

DDX 

~Jk 

3x 

GM 

body angle p 

DGP 

ar 

DGX 

il 

F 

body radius , f 

IMAX ' 


KMAX _ 

-see Description of Input 

Other 

program variables are 

FF 

2 

body radius squared, f 
2 

FFP 

3f 

3$ 

2 

FFX 

8f 

*2 2 

FFXP 

o r 
3x3<}> 
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FFXPP 


FFXXP 




Subroutine BGE0M 


Subroutine BGE0M reads the input data for the body geometry and calcu- 
lates the coefficients for the two-dimensional cubic spline function used 
to represent the body geometry. A call to subroutine BGE0M has the form 

CALL BGE0M ( IMAX , KMAX) 


where the arguments are defined in Description of Input. Other program 

variables are 
3 2 

AMD(I,K) 2-|- (x^) 

AMDD(l) ™1 temporary storage for elements of AMD(t,K) 
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AMF(I.K) 

AMFF(I) 

AMFF(K) 

AMFT(J) 

AMG(I.K) 

AMGG(I) ‘ 
AMGG (K) 

BMD(I,K) 

BMDD(I) 

BMF(I.K) 

BMFF(I) 

BMG(I.K) 

BMGG(I) 

DLP(I,K) 

DLPP(I) 
DLPP(K) . 
FN (I»K) 
FNN(I) 
FNK(K) _ 
FT(J) 

GAMU.K) 
GAMM(I) 
GAMH(K) _ 
HP(K) 
HPT(J) 


} 


3^^ 1 K 

temporary storage for elements of AMF(I,K) 

S 2 f 2 

— (40 corresponding to input coordinates 


3x 


-temporary storage for elements of AMG(I,K) 

3 2 

— 2 AMD(I.K) 

3x 


temporary storage for elements of BMD(I,K) 


3x 


2 AMF(I.K) 


temporary storage for elements of BMF(I.K) 
3 2 

— r AMGCI.K) 

rv 


temporary storage for elements of BMG(I,K) 
3 f2 , v 

JT ( W 

-temporary storage for elements of DLP(I f K) 

f2 <vV 

-temporary storage for elements of FN(I,K) 


E 2 

H~ (x i»v 


f corresponding to input coordinates 

2 


[-temporary storage for elements of GAM(1,K) 

*k ■ W 


- Vi 


HX(I) - x i _ 1 

I subscript corresponding to axial position 

J subscript corresponding to input circumferential position 

K subscript corresponding to rearranged circumferential position 

JM number of input circumferential locations 

JMM(l) number of input circumferential locations at x ■ X(I) 

PH(K) rearranged circumferential position, <)> k 

PHI ( J ) input circumferential position, <|>, 

X(I) input axial position, x. 

.2 

XMD;(I,K) • 2 -2 DLP (I ,K) 

Bx 

XMDD(I) temporary storage for elements of XMD(I,K) 

a 2 

XMF(I,K) — 5- FN(I,K) 
ax 

XMFF(I) temporary storage for elements of XMF(I,K) 

,2 

XMG(I,K) —y GAM(I »K) 

8x 

XMGG(I) temporary storage for elements of XMG(I,K) 

Y (J) input y-coordinate , see Description of Input 

YY(I,J) input y-coordinate at x = X(I) 

Z(J) input z-coordinate , see Description of Input 

ZZ(I,J) input z-coordinate at x = X(I) 

Note: The dimensions on the variables involving the subscripts I> J» and 

K in the program listing are 20 for each subscript. 

Subroutine SPLINE 

For an array of dependent variables, this subroutine solves the tri- 
diagonal matrix system of algebraic equations to determine the coefficients 
for the spline function. A call to SPLINE has the form 

CALL SPLINE (H,DE jK^AM^E ,G,N) 

where the dummy arguments are 
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H(J arra y ^presenting the difference between the J+l independent 

variable and the J independent variable 
DE(J) array of dependent variables 

K number of points used for the spline fit; J*1,2,...,K 

AM(J) array of coefficients for spline function 

E,G derivative (slope) of the dependent variable at J-l and J-K, respectively 
N indicator variable: 

(1) N - - 2 indicates AM(2)=AM(1) and AM(K-1)=AM(K) are used in 
place of specifying E and G. 

(2) N = - 1 indicates AM(1)«0 and AM(K)=0 are used in place of 
specifying E and G. 

(3) N = 0 indicates E is specified but AM(K-1)«AM(K) is used in 
place of specifying G. 

(4) N = 1 indicates both E and G are specified. 

Other program variables are 


Note: 


) array of constants in the system of algebraic equations 

/ 

subscript indicating position of dependent and independent variables; 
J-l, 2, .. . ,K 

^ array of parameters used in the algorithm for solving the system of 
algebraic equations 

This subroutine calculates the coefficients for spline functions in 
both the axial and circumferential directions. Therefore, the di- 
mension on the subscripted variables must be the number of points used 
in the axial or circumferential directions, whichever is greater. 

The computer program listed in Appendix A uses a dimension of 20. 


Function FSPL 


The function FSPL calculates the value of the dependent variable 
corresponding to a given value of the independent variable for any one of 
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the several spline functions. The function FSPL appears in the program in 
the form 

FSPL (AMI , AM2 ,X1 ,X2 ,Y1 , Y2 ,X) 

where the dummy arguments are 

AMI coefficient for spline function at XI 

AM2 coefficient for spline function at X2 

X value of independent variable, where XI <_ X < X2 

XI \ two successive members of the array of independent variables, used to 
X2 | determine the spline function, which bracket the current value of the 

J independent variable X, 

Yl dependent variable corresponding to XI 

Y2 dependent variable corresponding to X2 

Function DFSPL 

This function calculates the first derivative of the spline function 
computed in function FSPL. It appears in the program as 

DFSPL (AMI , AM2 ,X1 ,X2 , Yl , Y2 ,X) 

where the dummy arguments are the same as those in function FSPL. 

Function D2FPL 

This function calculates the second derivative of the spline function 
computed in function FSPL. It appears in the program as 

D2FPL (AMI , AM2 ,X1 ,X2 , Yl , Y2 ,X) 

and the dummy arguments are the same as those in function FSPL. 

Subroutine STGPR 

Subroutine STGPR computes stagnation properties for perfect gas or 
equilibrium air. A call to subroutine STGPR has the form 

CALL STGPR (PIN , TIN , VIN , RH0IN , EMIN , PS , RH0S ,HS,KG) 

where the arguments are 
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freestream pressure, p (lb/ft 2 ) 
freestream temperature, T^ (*R) 
freestream velocity, V w (ft/sec) 
freestream density, (slug/ft 3 ) 
freestream Mach number, M 

CO 

stagnation pressure, p (lb/ft 3 ) 

stagnation density, p g (slug/ft ) 

stagnation enthalpy, H (ft 2 /sec 2 ) 

s 

indicator variable; KG»0 for perfect gas or KG-1 for equilibrium air 
program variables are 
freestream enthalpy, h^ (ft 2 /sec 2 ) 
enthalpy aft of normal shock wave, h^ (ft 2 /sec 2 ) 
pressure aft of normal shock wave, p 2 (lb/ft 2 ) 
density aft of normal shock wave, p^ (slug/ft 3 ) 
gas constant for, ai r, R f (53. 35045 x 32.1741 ft 2 /sec 2 - e R) 
velocity aft of normal shock wave, v 2 (ft/sec) 

Subroutine PR0P 

This subroutine computes the local fluid properties at the edge of 
the boundary layer from the stagnation properties and local pressure 

for a perfect gas or equilibrium air. A call to subroutine PR0P has the 
form 

CALL PR0P (PS,HS,RH0S,PE,HE,RH0E,UE,VISC,EMS,KG) 
where the dummy arguments are 

rs 

PS stagnation pressure, p (lb/ft ) 

® 9 ‘ 9 

HS stagnation enthalpy, H g (ft /sec ) 

RH0S stagnation density, p g (slug/ft 3 ) 

PS local surface pressure, p (lb/ft 2 ) 

HE enthalpy at edge of boundary layer, h g (ft 2 /sec 2 ) 

RH0E density at edge of boundary layer, p (slug/ft 3 ) 

UE velocity at edge of boundary layer, U (ft/sec) 

Vise coefficient of viscosity at edge of boundary layer, (slug/ft-sec) 
EMS Mach number squared at edge of boundary layer, M 2 


PIN 

TIN 

VIN 

RH0IN 

EMIN 

PS 

KH0S 

HS 

KG 

Other 

HIN 

H2 

P2 

RH02 

RIN 

V2 
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KG indicator variable; KG-0 for perfect gas, KG-1 for equilibrium air 

Subroutine SPCH 

Subroutine SPCH computes the turbulent skin friction coefficient from 
the Spalding-Chi method based on momentum thickness. The corresponding 
turbulent heat transfer rate is calculated from von Karmen' s form of 
Reynolds Analogy. A call to subroutine SPCH has the form 

CALL SPCH ( ZW , HE , RH0E , UE , HS , RNM, CF , QW , QWB ,NCS ) 


where the arguments are 
ZW wall enthalpy ratio, 5 

w 2 2 
HE enthalpy at edge of boundary layer, h (ft /sec ) 

e 2 

RH0E density at edge of boundary layer, p (slug/ft ) 

UE velocity at edge of boundary layer, U (ft/sec) 

2 2 6 

HS stagnation enthalpy, H (ft /sec ) 

s 


momentum thickness Reynolds number, Rn^ 

turbulent skin friction coefficient, C- 

1 2 

turbulent heat-transfer rate, ^ (BTU/ft -sec) 

a „ ,/( h -h ), (ETU-sec/ft 4 ) 

htf.turb aw w 

indicator variable; if NCS-0 the value of U from the previous step 
is used in the Newton-Rhapson iteration scheme for U, if NCS«*1 the 
initial value of U for the iteration scheme is calculated from an 


approximate equation 


Other 

program variables are 


A0G 

h /h 
aw e 


A0S 

h /h 
aw w 


CFI 

incompressible skin friction coefficient, ^ 


FC 

F , see note below 
c 


FRD 



PR 

Prandtl number, Pr 


RAF 

Reynolds Analogy Factor, 

- 

REC 

recovery factor (0.89) 


S0G 

h /h 

* 


w e 
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u 


0.4 U+ 

Note: The equation used to compute FC in the program listing in Appendix 

A Is valid only for a perfect gas. After the listing of the program 
was printed, additions were made to subroutine SPGH to calculate FC 

for equilibrium air as indicated in Part I - Description of Basic 
Method. 

Subroutine PRESS 

Subroutine PRESS computes and prints modified Newtonian pressures 
corresponding to input body geometry coordinates, and then a two-dimensional 
spline fit to this data is obtained. If input pressure data are supplied, 
this pressure data is read-in and also printed in this subroutine, and a 
new two-dimensional spline fit to this data is calculated to replace the 

spline fit to the modified Newtonian pressures. A call to subroutine PRESS 
has the form 

CALL PRESS (KP.PNEWT) 

where the arguments are 
KP see Description of Input 


PNEWT (1 - Pa yp s ) 

Other program variables are 



3 2 p/p 


AMP(I,K) 

s 

3<J> 2 

<v* k > 


3 2 p/p 


AMPT(J) 

s 

3<j> 2 

<v 


BMP(I.K) ~ AMP (I ,K) 

3x 

BMPP(I) temporary storage for elements of BMp(l,K) 

FF body radius, f 

HX2 (I) X2(I) - X2(I-1) 

HPT(J) ^ 

H?2(K) $ k - ^_ 1 
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see Description of Input 

subscript corresponding to input circumferential position 

number of input circumferential positions 

number of input circumferential positions at X2(I) 

subscript corresponding to rearranged circumferential position 

see Description of Input 

indicator variable; LL-1 for modified Newtonian pressures at 
input body geometry coordinates, and LL«2 for input pressure data 
input circumferential position, 
rearranged circumferential position, ^ 


PN(I,K) p/p 8 (* ± ,+ k ) 


PNN(I) 

PNN(K) 

P0P 

PR 

PT (J) 


-temporary storage for elements of PN(l,K) 

pressure ratio, p/p g 
pressure ratio, p/p 

s 


PX(K) 


(0, PH2 (K) ) 


PXE(K) 


(X2 ( IM2X) , PH2 (K) ) 


PXPP(K) 


PX(K) 


PXPPE(K) 


PXE(K) 


XMP(I,K) 


a 2 p/p c 


(x i*V 


XMPP(l) temporary storage for elements of XMP(I,K) 

XT axial location, x 

X2(l) /XT 

Y(J) y-coordinate of 

YY(I,J) y-coordinate at x^ , <j> ^ 
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Z(J) z-coordinate of $ 

ZZ(I,J) z-coordinate at 

J 

Subroutine PDER 

For a given position on the body, this subroutine calculates the 
pressure ratio (p/p s > and its derivatives from the two-dimensional spline 
pressure function. If pressure data are input this function is a spline 
fit to the input pressure data; otherwise, the function is a spline fit to 
the modified Newtonian pressures corresponding to the input body geometry 
coordinates. A call to subroutine PDER has the form 

CALL PDER (XX , PPH , PE , DPDP ,DP2P , DPDX , DPXP , DPXX) 

where the dummy arg laments are 

XX axial position on body, x 

PPH circumferential position on body, $ 

PE pressure ratio, p/p 

9(P/P 8 > 

T\TJT\Ti O 


3 .( P/P 8 > 


3 (P/Pg) 


3 (P/P ) 

o 


3 (P/P s > 


Other program variables are defined in subroutine PRESS 


Subroutine LAGHN 


This subroutine calculates the axial pressure derivatives (}with re- 
spect to /x. ) at the nose (x«0) and end of the body for <j> » <f>, from 
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Lagrangian Interpolation functions. These pressure derivatives are used 
as "end conditions" for pressure spline functions in the axial direction in 
subroutine PRESS. A call to subroutine LAGRN has the form 

CALL LAGRN (XX,X,FN,IMAX,M,N ,F,FP) 


where the dummy arguments are 

XX axial location where derivative is to be evaluated, vx 

X(I) array of independent variables (v^x) 


XX 

X(I) 

FN(J) 

IMAX 

M 


array of dependent variables (p/p g ) 
number of axial stations 

points to the left of XX to be used in Lagrangian Interpolation 
points to the right of XX to be used in Lagrangian Interpolation 
p/p at x«XX and 


3(P/P S ) 


at x=XX and 

Subroutine TRANS 


Subroutine TRANS is used only if transition is specified by geometric 
location on the body (KTR«-1 in input data) . It reads the body coordinates 
specifying the beginning and end of transition, and it calculates the co- 
efficients for a one-dimensional spline function for each set of data. A 
call to subroutine TRANS has the form 

IF(KTR .EQ. -1) CALL TRANS 


and the 

program variables are 
32 

AMTB(L) 

9<j) 2 

9 

XB(L) 

AMTE(M) 

3 

3<j> 2 

XE(M) 


subscript corresponding to a body coordinate used to specify 
beginning of transition; L=l,2, . . . ,LMAX 

number of coordinates used to specify beginning of transition 
(see Description of Input) 
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subscript corresponding to a body coordinate used to specify end 
of transition; M«l,2,. . . ,MMAX 

number of coordinates used to specify end of transition (see 
Description of Input) 

circumferential location of input coordinate for beginning of 
transition, 4 . <L> 

circumferential location of input coordinate for end of tranr 
sition, 4 (M) 

\ 

Cartesian coordinates of position L specifying beginning of 
| transition 

^ Cartesian coordinates of position M specifying end of transition 

Subroutine BEGTR 

This subroutine is used only when transition is specified by geometric 
location on the body (KTR=-1 in input data) . For a given circumferential 
location (PPH) it calculates the corresponding axial location (XTB) of the 
beginning of transition. A call to subroutine BEGTR has the form 

IF(KTR .EQ. -1) CALL BEGTR (PPH, XTB) 

Subroutine ENDTR 

This subroutine also is used only when transition is specified by 
geometric location on the body (KTR = — 1 in input data) . For a given cir- 
cumferential location (PPH) it calculates the corresponding axial location 
(XTE) of the end of transition. A call to subroutine ENDTR has the form 

IF(KTR ,EQ. -1) CALL ENDTR(PPH.XTE) 

DESCRIPTION OF INPUT 

The input data for the computer program are listed below in the proper 
sequence. Read statements for the input data appear in the Main Program 


M 

MMAX 

PHB(L) 

PHE (M) 

XB(L) 

YB(L) 

ZB(L) 

XE(M) 

YE(M) 

ZE(M) 
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unless noted otherwise. In a subsequent section the input data used in the 
examples in Part I of this report are listed. The steps which describe the 
input are: 

Step I. 

READ (1,1000) NAM 
1000 F0RMAT (20A3) 

NAM is the name used for the body designation. 

Step 2o 

READ (1,1) IMAX, KMAX 
1 FORMAT (215) 

This read statement appears in subroutine BGE0M. 

IMAX number of axial stations used to input body geometry data, including 
the nose point (x“Q) and the last axial station. Although IMAX in- 
cludes the nose point (x=*0) , body geometry data are not input for 
this point since the body radius f is zero for all circumferential 
angles (See Sketch 1 on next page) 

KMAX number of equally spaced circumferential planes to be used to spline 
fit the body geometry in the circumferential direction, including 
the windward (<j>«0) and leeward (<j>“TT) planes but only for the half 
body 0<j><7To These equally spaced circumferential planes are gener- 
ally different from the circumferential angles corresponding to the 
input coordinates described in step 3. Although the input body co- 
ordinates are not generally equally spaced in the circumferential 
direction and generally differ from one body cross section (axial 
station) to another, a one-dimensional spline fit to the input data 
at each cross section is used in the program to calculate rearranged 
body points which are equally spaced in the circumferential direction. 
This rearrangement of geometry data is needed to facilitate the spline 
fit in the axial direction. (See Sketch 1.) 

Dimension statements in the program require that 
IMAX £ 20 and KMAX £ 20. 

) 
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Sketch 1 - Body geometry 

Step 3 . / 

D0 11 1=2 , IMAX 
KEAD(1,3)JM,X(I) 

3 F0|WAT(I5,F1O.5) 

BEAD(1,4) ((Z(J) ,Y(J)) , J=l, JM) 

4 F0RMAT(8F1O»5) 

11 C0NTINUE 


K=MAX 



\ 


y 

View from Rear of 
Rearranged Circumferential 
Planes (Equally Spaced) 


The statements above appear in subroutine BGE0M, and they are the 
read statements to input the Cartesian coordinates of the points used to 
describe the body geometry . 

I subscript corresponding to axial station number, where 1=1 is the 
nose point (x=0) and I = IMAX is the last axial station number. As 
noted in step 2, body geometry data are not input for the nose 
point (1=1) » and therefore the D0 loop listed above starts with 1=2. 
JM number of circumferential points input at axial station number I 
(see Sketch 2 on next page.) 

X(I) axial distance from the nose point (x-0) to axial station I. The 

axial stations are not required to be equally spaced, but the spline 
fit in the axial direction is generally more accurate when the axial 
stations are spaced approximately equally® (See Sketch 1 above.) 


2 , 


J subscript corresponding to an input circumferential point, where 
J=1 is the windward plane and J=JM is the leeward plane. 

Z(J) z-coordinate for point J 
Y(J) y-coordinate for point J 


J»JM 



y 

Sketch 2. - Input coordinates of body cross section (viewed from rear). 

The points J=l,2,3, . . . ,JM do not have to be equally spaced. However, it 
is required that J=1 corresponds to the windward plane, J=JM corresponds 
to the leeward plane, and 4>(J+l)><j> (J) where <)>(J)=tan (Z(J)/Y(J)). Also, 
JM may vary from one axial station to another and is not related to KMAX 
except the dimension statement in the program has JM<_20. The coordinates 
X(I), Y(J) , and Z(J) may be input with any consistent unit of length since 
the constant CST in step 4 below is used to convert this unit to feet. 

Step 4. 

READ (1,900) PIN, TIN, VIN,TRI , THE, CST ,HMAX,KG,KTR,KS ,NMAX,KP0, 

900 F0RMAT(7E1O. 0/5110) 

2 

PIN freestream static pressure, p^ (Ib/ft ) 

TIN freestream temperature, T^ (°R) 

VIN freestream velocity, (ft/sec) 
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TRI value of integrated unit Reynolds number Rn^ or momentum thickness 
Reynolds number Rn ffi , if either of these two parameters is used to 
specify the beginning of transition. If transition is specified by 
geometric location or not specified at all, any number may be input 
for TRI 

TRE same explanation as given above for TRI except TRE corresponds to 
the end of transition 

CST constant to convert units of input body geometry coordinates to feet 
HMAX maximum step size along a streamline for integration by Runge-Kutta 

method; HMAX has the same units as the input body geometry coordinates 
KG indicator variable; KG=0 for perfect gas, KG-1 for equilibrium air 

KTR indicator variable for transition; KTR=-1 for beginning and end of 

transition based on geometric location, KTR-0 for transition based 
on integrated unit Reynolds number Rn^ or KTR-1 for transition based 
on momentum thickness Reynolds number Rn^. When transition is not 
specified, KTR may be any other integer as long as KS=0 is input 
KS KS=0 for laminar heating rates only, otherwise any other integer may 
be input 

NMAX maximum number of integration steps along a streamline 
KP0 when the additional print out (see Description of Output) is desired, 
set KR0«=1 ; otherwise any other integer may be input 
Step 5. >' 

The data in this step are input only if the beginning and end of 
transition is specified by geometric location (KTR=-1). Otherwise, skip 
this step and go to step 6. 

The statements below appear in subroutine TRANS, and they are the 
read statements for the Cartesian coordinates of points along the body used 
to specify the beginning and end of transition by geometric location. 
READ(1,1)LMAX 
1 F0RMAT(I5) 

READ(1,3) ((XB(L) ,yb (l) ,ZB(L)) ,L*1,LMAX) 

3 FORMAT (8E10.0) 

RE|JD(1,1)MMAX 

e 

READ (1 , 3) ( (XE (M) ,YE (M) , ZE (M) ) ,M-1 ,MMAX) 
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The variables are 

LMAX total number of input points used to specify beginning of transition, 
including one in the windward plane and one in the leeward plane 

(0 <$<].}) 

L subscript corresponding to input point to specify beginning of tran- 
sition; where 1<L<1MAX. , L«1 must be the windward plane (((>“0) and 
L«LMAX must be the leeward plane ((j>“ 7 r) 

XB(L)J x, y , and z coordinates, respectively, for input points used to 
YB(L)\ specify beginning of transition by geometric location* The units 
ZB(L) are the same as those used for specifying the body geometry coordi- 
nates. Also it is required that (L-hl) > <f> (L*) where 
<KL)«tan -1 (ZB(L)/YB(L)) and (t>(l)»0, <KLMAX)=tt 
T he variables listed below pertain to the end of transition, and the same 
general comments given above for the variables pertaining to the beginning 
of transition apply to these variables also* 

MMAX total number of input points used to specify end of transition 
M subscript corresponding to input point to specify end of transition, 
1<M£MMAX 

x, y, and z coordinates, respectively, for input points used to 
specify end of transition 

Sdep 

RE AD ( 1 , 2 ) ALP D , ZW , KBM , KP 
2 FORMAT (2F10 o 5,215) 

ALPD angle of attack, a (degrees) 

ZW ratio of wall to freestream stagnation enthalpy, £ w 

KBM total number of streamlines to compute for the half body (0<jj><jr) 

KP indicator variable; KP«0 for simplified streamlines (see Part I for 
explanation) , KP=X for streamlines computed from modified Newtonian 
pressures, KP-2 for streamlines computed from input pressure distri- 
bution 
Step 7 0 

The data in this step are input only if surface pressures from some 
other source are supplied (KP 5 ^) . Otherwise, skip this step and go to step 


XE (M) 
YE (M) 

ze(m) 
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8. The pressure data are read-in as the ratio of surface pressure to stag- 
nation pressure around the circumference of the body at several axial 
stations, in a manner similar to the body geometry coordinates described 
above. The statements below appear in subroutine PRESS. 

READ(1,15)PN(1,1) ,IM2X,KM2X 
15 F0RMAT(FlO.5i2I5) 

D0.11 I-2.IM2X 
IF(LL„NE. 1)READ(1,3) JM,XT 

3 F0RMAT(I5,F1O.5) 

IF(LL.NE.1)READ(1,4) ((Z(J) ,Y(J) ,PT(J) ,J=1,JM) 

4 F0RMAT(6F1O.5) 

11 C0NTINUE 


IM2X 


KM2X 


The program variables are: 

PN(1,1) ratio of pressure at the nose point (x**0) to stagnation pressure 
number of axial stations where pressure data are supplied, in- 
cluding the nose point (x“0) and the last station. Note that 
these axial stations are not required to coincide with the axial 
stations used above for the input body geometry coordinates 
number of equally spaced circumferential planes to be used to 
spline fit the rearranged pressure data in the circumferential 
direction, including the windward (<|>=0) and leeward planesT 

Although the value of KM2X generally differs from KMAX, the same 
comments given for KMAX above apply here also, 
axial station number, where I«*l is the nose point (x«0) and 
I=IM2X is the last station number 

number of circumferential positions where pressure data are input 
axial distance from the nose to the local axial station, in the 
same units as input body geometry coordinates 

y and z-coordinates , respectively, of position J where pressure 


JM 

XT 


Z(J) 

Y(J) 


PT(J) 

J 


ratio is input (units are same as input body geometry coordi- 
nates). Note that $(J+1)><KJ) , where <j>(J)=tan“ 1 (2(j)/Y(J)) 
ratio of surface to stagnation pressure at position J 
subscript corresponding to a circumferential position where pres- 
sure ratio is input, where J-l is the windward ($=0) plane and 
J=JM is the leeward (<)>=* 7 r) plane. 
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Note that JM may vary from one axial station to another and the points 
J=l,2, e * c ,JM do not have to be equally spaced* However, the dimension 
statement in the program limit JM, IM2X and KM2X to 20* 

Step 8 0 

KB=1 

g 

3 RE AD (1,4) BET AD 

4 F0RMAT(F1O.5) 

19 IF (KB * GE * KBM) G0 T0 10 
KB«KB4*1 
G0 T0 3 

The statements above read-in a value of BETAD for each streamline to be 
computed, where 

KB streamline number, 1<KB<KBM 

KBM total number of streamlines for half-body 0 <$<tt 

BETAD angular orientation (2) of each streamline on the e-circle (see 

page 36 of Part I), measured in degrees with BETAD=0. corresponding 
to windward streamline, BETAD “ 180. corresponding to leeward stream- 
line, and Oe<BETAD<180, A total of KBM values of BETAD must be input. 

DESCRIPTION OF OUTPUT 

"Ah example of the output obtained from the computer program for two 
streamlines (g^l 0 and 45^) on a blunted 15’ half-angle cone at 20° inci- 
dence and » 10*6 is given in Appendix C. The format for the output was 
changed slightly in Appendix C to facilitate printing it in this report. 

The output data printed are: 

1. Body designation specified by NAM in the input data 

2c If transition is specified by geometric location (KTR^-l) , the 

Cartesian coordinates of the input points for the beginning of tran- 
sition (XB,YB,ZB) and end of transition (XE,YE,ZE) are printed. The 
total number of points used to specify the beginning of transition 
is given by LMAX, and the number for the end of transition is MMAX. 

3. Freestream properties: pressure (p w ) » temperature (T^) , density 

(p ) , velocity (V ), and Mach number (M )* 

4* Stagnation properties: pressure (p ), density (p ), and enthalpy (H ). 

s s s 
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5. 

6 . 

7. 

8 . 


10. 


Input parameters TRI,TRE,CST,KG,KTR,KR,HMAX,NMAX,KP0,KS,IMAX, and KMAX. 
The input body geometry coordinates are printed, and the modified 
Newtonian pressure (ratioed to stagnation pressure) is printed next 
to the coordinates of each point. The modified Newtonian pressure 
ratio corresponding to each input body coordinate is printed even 
when input pressure data from some other source is supplied. 

If input pressures (ratioed to stagnation pressure) are supplied, this 
data is printed along with the corresponding body coordinates. 

ANG. OF ATTACK - angle of attack, ct(degrees) 

ZETA(WALL) - wall enthalpy ratio, £ _ 

NO. OF STREAMLINES TO BE CALCULATE - KBM 

B - ratio of principal radii of curvature at stagnation point, R /R n 

TW - wall temperature, T (°R) 

w 

QWS - stagnation-point heat-transfer rate, q (BTU/ft^-sec) 

HSP - QWS (6006) /(H g -h w ), (BTU/'f t 2 -sec-*R) whlre the factor 6006 is 

the specific heat at constant pressure for a perfect gas (air) , 

Hg is the stagnation enthalpy, and h w is the wall enthalpy. 

Note that for a perfect gas HSP is the stagnation-point heat- 
transfer coefficient. 


For each position along a streamline, the normal print-out gives the, 
following data: 

X - x-location from nose | 

Y y-location I same iin-ft*® ae 


Z - z-location 


same units as 
input body 
geometry data 


,, '-io. L.O 

H - scale factor or ’ equivalent radius 11 , h 
S - distance along the streamline, measured 
from the stagnation point J 

BBAR - pressure gradient parameter, g" 

P/PS - ratio of local to stagnation pressure, p/p 

QW/QWS - ratio of local to stagnation-point heat-transfer rate, q^/q^ 

If KP0«1 is input, then the following additional data are printed at 
each position along a streamline: 


UE - velocity at edge of boundary layer, (ft/sec) 

HE - enthalpy at edge of boundary layer, h g (ft 2 /sec 2 ) 
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3 

RH0E - density at edge of boundary layer, p (slug/ft ) 
MUE - coefficient of viscosity at edge of boundary layer, 
(slug/ft-sec) 

MACH NO - Mach number at edge of boundary layer, M 

e 

RN/FT - Reynolds number per foot, pH /y (ft ) 

e e e 

MOM THK - momentum thickness, 6 (ft) 

in 

RNM - momentum thickness Reynolds number, Rn 

m 

RNl - integrated unit Reynolds number, Rn^ 

R - body radius, f (in units of input body data) 

PHI - angular position on body, <J> (degrees) 


DR/DPHI - 

D2R/DPHI2 - 

DR/DX - 
D2R/DX2 - 


3 f 

*rr (same units as R) 
0<J) 

3 2 f 

— — (same units as R) 

3<j> 


3f 

3x 

3 2 f 

3x 2 


(units of R 


QB/QSB - (QW/QWS)*(H s -h w )/(h aw -h w ) where H g is frees tream stagnation 
enthalpy, h is adiabatic wall enthalpy, and h^ is wall 
enthalpy o 


DIAGNOSTIC MESSAGES 

Listed below are descriptions of diagnostic messages printed in the 
output : 

1 8 NEGATIVE ANGLE OF ATTACK ENCOUNTERED, STAG. POINT SET TO 0.0, PROGRAM 
COUNTINUING. If a negative angle of attack a is input to the program, 
then the program places the stagnation point at the nose (x=0) in 
subroutine STAGN , assigns the value -1 to the indicator variable ISTG, 
prints the message above and continues the computations* For those 
examples where a negative angle of attack is desired, the body geome- 
try data should be input upside down and then a positive angle of 
attack can be used in the input* 
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2. STAGNATION POINT NOT FOUND, RUN ABORTED. An iterative technique is 
used to locate the Newtonian stagnation point in subroutine STAGN. 

If the stagnation point should not be located, the program assigns 
ISTG-0, prints the message above, and terminates the calculations. 

3. CALCULATIONS FOR BETA - HAS TERMINATED AFTER INTEGRATIONS. 

If the computations along a particular streamline should reach the 
end of the body before the number of integration steps (L) has reached 
the maximum allowed (NMAX in the list of input) , then the message 
above is printed at the end of the computations for that streamline. 

4. *** SUBSCRIPT OUT OF RANGE ***, CALCULATIONS FOR STREAMLINE NO. 
DISCONTINUED. When the number of integration steps (L) along a 
streamline has reached the maximum allowed (NMAX) before the end of 
the body is reached, the message above is printed. 

5. INTEGRATION OF STREAMLINE DISCONTINUED DUE TO STEP SIZE, 

*** H = ***» When the integration step size (H) in function 

RUNGE is decreased to a value less than 10~ 8 , the message above is 
printed and computations for this streamline are discontinued. 

6. SCALE FACTOR (H) GOING NEG. , *** H = STREAMLINE DISCONTINUED. 

The message above is printed if the scale factor or "equivalent 
radius (H) becomes less than 10 and the computations for that 
streamline are discontinued. 

7. REY NO (RNM) GOING NEG, *** RNM *= , STREAMLINE DISCONTINUED. 

If the momentum thickness Reynolds number (Rn ) becomes less than 
—4 m 

10 , the message above is printed, and the computations for that 

streamline are discontinued * 

8. **** BEGINNING OF TRANSITION ****. This message is written just before 
the print— out at the beginning of the transition region. 

9. **** END OF TRANSITION ****,, This message is written just after the 
print-out at the end of the transition region. 

TEST CASES 

Test cases were calculated for: Case I. - A blunted 15‘ half-angle 

cone at ot - 20 e and - 10.6, Case II. - A blunted 70° slab delta wing 
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at a - 10 e and M - 8, and Case III. - The HL-10 Lifting Body at a - 20* 

00 

and M -10. Inputs for these cases are given below, and output data for 

00 

two streamlines on the blunted 15° half-angle cone are given in Appendix C. 
Results for all three cases are presented graphically in Part I of this 
report | and the results for the blunted 15 c half-angle cone are also shown 
to compare favorably with experimental data. 

Case I. - Blunted 15® Half-Angle Cone 

Input data for a blunted 15° half-angle cone at a •» 20® and « 10.6 
are presented on the pages which follow. Written beside each section of 
input data is the step number which describes that data in Description of 
Input o Step 1 corresponds to the body designation, and for Step 2 the body 
geometry was specified by input coordinates at 19 axial stations (IMAX»19) 
and 20 planes around the half body for the rearranged circumferential points 
(KMAX«20) . As noted in Part I of this report, the body geometry could be 
represented by simple analytical expressions, but the coordinates of 20 
points around the circumference of the half body at 19 axial stations were 
used to generate the body geometry in order to test the accuracy of the 
doubly cubic spline function. This data corresponds to Step 3 of the in- 
put, and the coordinates are input in inches. 

2 

Frees tream properties are shown in Step 4, where p^ “ 2.6614 Ib/ft » 

- 89.971 °R, and V w - 4928.1 ft/sec. Since only laminar heating rates 
were calculated here, TRI=0 and TRE=0 * In order to convert the Input 
coordinates from inches to feet, a value of CST-0. 0833333 was used. The 
maximum step size was set at HMAX»0 o 4, and a perfect gas (I<G®0) was 
specified. Again, since only laminar heating rates were calculated, KTR 
could be any integer (KTR“0 was used here) and KS=0. The maximum number 
cf integration steps along a streamline was input as NMAX=*100, and since 
the additional print-out along each streamline was desired, KP0=1. 

Step 5 of the input data was not included because transition was not 
required. For step 6, an angle of attack (ALPD) of 20°, Wall enthalpy 
ratio (ZW) of 0.251, number of streamlines (KBM) of 2, and streamline com- 
putation by modified Newtonian pressures (KP^l) were input. Step 7 of 
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the input data was omitted because modified Newtonian pressures were used 
in lieu of an input pressure distribution* 

The input data given here only have two streamlines to be calculated 
for Step 8 - BETAD of 1 degree and 45 degrees. Note, however, that the 
data presented in Part I of this report used BETAD values of 1 # , 10*, 15°, 
and 45* . 

A complete listing of the output data for this case is given in Appen- 
dix C. For an explanation of this data, see Description of Output. From 
the data listed in Appendix C, it is observed that the integration step 
size along a streamline is small near the stagnation point, but it in- 
creases rapidly downstream due to the variable step-size Runge— Kutta 
integration routine. The streamlines move from the windward! side of the 
body and approach the leeward plane downstream. The scale factor or 
"equivalent radius" increases all along the BETAD “ 1° streamline, but 
it decreases along the BETAD « 45° streamline after it has "essentially" 
reached the leeward plane (PHI ■ 180.). It should be noted, however, that 
the streamlines are continued into the "shadowed" region of the body where 
the surface pressure is assumed to be freestream static pressure. It is 
possible that the flow separates from the body somewhere in thiB region, 
and therefore the calculated properties are questionable there. Along 
each streamline the heating-rate ratio QW/QWS decreases and them gradually 
levels out. The computational time (CPU) for this case on an IBM 360/75 
computer was 1:47 minutes. 

Case II. - Blunted 70* Slab Delta Wing 

Input data are given on the following pages for a 70° swept delta wing 
with a cylindrical leading edge which is tangent to a flat slab on the 
upper and lower surfaces. The blunt nose is a spherical cap, and the 
radius of both the spherical cap and cylindrical leading edge is one foot. 
Again, the step number corresponding to that in the Description of Input 
is written beside each section of the input data. 

Step 1 corresponds to the body designation, and in Step 2 the body 
geometry is specified by input coordinates at 19 axial stations (IMAX“19) 
and 19 planes around the half body for the rearranged circumferential points 
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Input for Case II (concluded) 







(KMAX-19). The coordinates of 13 points around the circumference of the 
half body at 19 axial stations were used to generate the body geometry in 
Step 3. All body coordinates were input in feet* 

3 2 

In Step 4 the freestream properties are p “10 lb/ft , T « 416°R, 

and V ■ 8000 ft/sec« For this case transition was assumed to begin at a 

00 3 3 3 

momentum thickness Reynolds number of 10 (TRI*10 ) and end at 2 x 10 

3 

(TRE“2 x 10 ) for each inviscid surface streamline. Since the input body 
coordinates have units of feet, CST -1.0. A maximum integration step 
size of HMAX * 0*05 was input, and a perfect gas (KOO) was specified. 

Since transition was based on momentum thickness Reynolds number, KTR-0 
and KS«1. The maximum number of integration steps allowed along a stream- 
line was NMAX«*100, and KP0-1 was input to obtain the additional print-out 
data along each streamline. 

Input data corresponding to Step 5 were omitted because transition was 
not specified by geometric location. In Step 6 the angle of attack was 
ALPD*10°, the wall enthalpy ratio was set at ZW«0.4, ; the number of stream- 
lines to be calculated was KBM-8, and the streamlines were calculated from 
modified Newtonian pressures (KP“1) . 

Step 7 was omitted since surface pressures from some other source were 
not used. For the 8 streamlines to be calculated, the values of BETAD in 
Step 8 are 0, 1, 10, 30, 60, 90, 170, and 179 degrees. 

The computational time (CPU) for this case was 3:02 minutes on an IBM 
360/75 computer* Surface streamlines and heating rate ratios are presented 
graphically in Part I of this report. 

Case III. - HL-10 Lifting Body . 

The third, and last, test case computed was for the HL-10 lifting body, 
which Is described in Part I of this report. Input data for this case are 
listed on the following pages, where again the Step numbers corresponding 
to the steps in the Description of Input are written beside each section 
of input data. Step 1 gives the body designation, and Step 2 lists the 
number of axial body stations (IMAX-20) and rearranged circumferential 
planes (KMAX«20) used to describe the body geometry. In Step 3 the 
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Input for Case III 
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Input for Case III (continued) 
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Input for Case III Concluded; 




coordinates of the input body points are listed at 20 axial locations and 
19 circumferential positions, where all coordinates are input in inches. 


The freestream properties are given in Step 4, where p 
416 R, and V_ »« 10 ft/sec. In this case transition was specified 


10 3 lb/ft 2 . 


by geometric location; therefore TRI and TRE could be any numbers (values 
of TRI-10. and TRE=20.were arbitrarily input). Since the input body co- 
ordinates were in inches, a value of GST - 0.08333 was used. The maximum 
integration step size was choosen to be HMAX « 18.1235, which is one-tenth 
of the distance from the nose to the last axial body station. Equilibrium 
air (KOI) was used for this case, and KTR=-1 was input because transition 
is specified by geometric location. For laminar, transitional, and turbu- 
lent heating, any non-zero value of KS may be used, and here KS-1 was input. 
The number of integration steps along a streamline was limited to NMAX“100, 
and only the normal print-out along each streamline was desired (KP0=O) . 

In Step 5 data was input to specify beginning and end of transition 


by geometric location. Three points were used to specify the beginning 
(IMAX-3) and end (MMAX-3) of transition. In each case the axial location 
was arbitrarily chosen to be x = 33.073" for beginning and x = 54.239" for 
end of transition. 

For Step 6, the angle of attack was ALPD - 20”, wall enthalpy ratio 
ZW ■ 0.1, number of streamlines KBM - 5, and the streamlines were calcu- 
lated by the "simplified method" discussed in Part I (KP=0) . Step 7 was 
omitted. For the five streamlines to be calculated, values of BETAD in 
Step 8 are 0, 45, 90, 135, and 180 degrees. • 

Heating-rate ratios for this case are presented graphically in Part I 
of this report for the windward streamline. The computational time for 
this case (CPU) was 1:48 minutes on the IBM 360/75 computer. 


ACCURACY AND LIMITATIONS OF METHOD 

The method used in the computer program described herein is applicable 
to a large class of configurations subject to the following restrictions: 

1. The body must have at least one plane of symmetry, and the free- 
stream velocity vector must be parallel to this plane of symmetry. 
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2. The body radius r * f(x,<j>) must be single valved. 

3* The nose of the body must be blunted, and at x«0 and 

A two-dimensional cubic spline function is utilized to represent both 
the body geometry and the input pressure distribution. The computer program 
listed in Appendix A has dimension statements which limit the number of 
axial and circumferential positions to 20; however, these dimensions could 
be changed to accommodate as many positions as desired. As indicated in 
Part I of this report, the spline function does not smooth out data points 
(although it does smooth derivatives), and therefore data points used in the 
spline function should be reasonably accurate. When pressure data from some 
other source are used, they should be smoothed before using them in the com- 
puter program. 

The accuracy of the spline function is affected by both the number and 
location of data points used. More accurate results are generally obtained 
when the data points are nearly evenly spaced in the axial and circumferen- 
tial directions. In some cases a large number of data points are needed to 
make the spacing between data points nearly even and also have a reasonable 
number of points in regions where the data is highly varying. As an example, 
consider the cross sectional geometry of a blunted slab delta wing far down- 
stream of the nose. At this axial station the thickness of the wing is small 
compared to the span, and a large number of body coordinates are needed to 
make the data evenly spaced in the circumferential direction and still have 
a reasonable number of points near the tip of the wing. 

Input freestream properties for this method are limited to supersonic 

speeds, and when modified Newtonian pressures are used the speed should be 

hypersonic. Although the angle of attack is not limited to small angles, 

heating rates on the lee side of the body are questionable at large angles 

of attack. The approximate methods used to calculate heating rates are 

more accurate for cold walls (£,<<1) than hot walls. 

w 

In Part I of this report the present method was found to predict 
laminar heating rates very well on the blunted 15° half-angle cone at 20 fc 
incidence and « 10.6. Additional comparisons of the present theory with 
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experimental data on other body shapes and for transitional and turbulent 
heating are needed to assess the accuracy of the theory more thoroughly. 

CONCLUDING REMARKS 

The computer program described in this report was developed to calcu- 
late inviscid surface streamlines, and laminar, transitional, and turbulent 
heating rates on shuttle-type configurations. Approximate methods were 
employed to keep the computational time and storage down to relatively 
small values but provide reasonably accurate results. Input and output 
data were programed to be easily used by the ordinary engineer who is not 
trained as a computer specialist. Although the computer program was devel- 
oped on an IBM 360/75 computer, it is also compatible with the IBM 370/165, 
CDC 6400, and CDC 6600 computers. 
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APPENDIX A 

PROGRAM LISTING 
from IBM 360/75 Computer 

Contents 


Main Program - - 
Function RUNGE - 
Subroutine STAGN 
Subroutine NEWP 
Subroutine STRMLN 
Subroutine BGE0M 
Subroutine SPLINE 
Function FSPL 
Function DFSPL - 
Function D2FPL - 
Subroutine STGPR 
Subroutine PR0P 
Subroutine SPCH 
Subroutine PRESS 
Subroutine PDER 
Subroutine LAGRN 
Subroutine TRANS 
Subroutine BEGTR 
Subroutine ENDTR - 
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APPENDIX B 

DETAILED FLOW DIAGRAM 

Symbols used in the flow diagram conform to the International Organization 
for Standarization (ISO) Recommendation on Flowchart Symbols for Information 
Processing, and are consistent with the fewer symbols adopted by the American 
National Standards Institute, Inc. (ANSI). 

CONTENTS 


Main Program - - 
Function RUNGE - 
Subroutine STAGN 
Subroutine NEWP - 
Subroutine STRMLN 
Subroutine BGE0M 
Subroutine SPLINE 
Function FSPL - - 
Function DFSPL - 
Function D2FPL - 
Subroutine STGPR 
Subroutine PR0P - 
Subroutine SPCH - 
Subroutine PRESS 
Subroutine PDER - 
Subroutine LAGRN 
Subroutine TRANS 
Subroutine BEGTR 
Subroutine ENDTR 
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—0^19633 

P/PSC 151= 

;l - 

Z C 161= 

0.17804 

YC 161 = 

-0.22875 

P/PSC 161= 


ZC 171= 

0.13796 

YC 171= 

-0.25494 

P/PSC 171= 

j 

ii 

1! 

ZC 181= 

0.09412 

YC 181= 

-0.27417 

P/PSC 181= 


0- 76857 
0.71955 
0.66862 
0.61757 
0c 56799 
0.52106 
0.47797 
0. 43963 
0.40667 
0.37929 
0.35798 
0.34271 
0.33354 
0.33046 


0.74217 
0.73604 
0.71786 
0.68867 
0.65001 
0.60356 
0.55163 
0. 49664 
0.44069 
0.38611 
0.33439 
0.28704 
0.24517 
0.20924 
0.17952 
0.15593 
0.13810 
0.12572 
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2 1,19 1 = 

Q. 04771 

¥( 19)= -0.28592 

P/PS? 191 = 

0.11851 

78201= 

O.G 

¥E2Gi= -0.28987 

P/PSf 201 = 

0.11611 

:«00- NEW? . PRESS. AT 

INPUT COORD. 




•; FOR UM = 

20 AND X( 

41= 0. 

20564 



?**E Z,Y 

AND P COORDINATES ARE: 



,11 = 

0.0 

¥8 11= 

0.33458 

P/PSf 15= 

0.53630 

It 21= 

0^05507 

Y( 25 = 

0.33002 

P/PSf 25= 

0.53028 

•If .31 = 

0.10864 

¥1 .35 = 

0.31645 

P/PS8 .35= 

0.51260 

21 41= 

0.15924 

¥8 41= 

0.29425 

P/PSf 45= 

0.48416 

. «. .-55* 

0.20550 

Y8 51= 

0.26403 

P/PSf 55= 

0.44691 

Z! :6|= - 

0.24616 

YC 65= 

0.22660 

P/PSI 65= 

0.40308 

2! 7J= 

0.28010 

Yl 71= 

0.18300 

P/PS8 75= 

0.35501 

Z8 81= 

0.30640 

YC 85 = 

0.13440 

P/PSf 85= 

0.30460 

2( 91= 

0.32434 

Y8 95 = 

0^08213 

P/PSf 95= 

0.25510 

28105 = 

0.33344 

¥1105= 

0.02763 

P/PS8 105= 

0.20810 

ZflI5 = 

0.33344 

Y{ 115= 

-0.02763 

P/PSC1I1= 

0.16553 

2! 121= 

0.32434 

Yi 125= 

-0.08213 

P/PS< 121= 

0.12866 

71131= 

0.30640 

YC 135= 

—0^13440 : 

P/PSf 135= 

0.09776 

ii 21141= 

0.28010 

Y 8.141= 

-0.18300 

P/PSf 141= 

0.07320 

£5151= 

0.24616 

YC 151= 

-o.; 22660 

P/PSf 15) = 

0.05427 

* 211-65= 

0.20550 

¥8165= 

-0.26403 

P/PSI 165= 

0^04057 

21175= 

0.15924 

¥8175= 

-0.29425 

P/PSf 171= 

0.031 27 

2818}= 

0.10864 

¥1181= 

—0.31665 

P/PS 8185= 

0.02543 

28191= 

0.05507 

¥(191= 

-0.33002 

P/PS819I* 

0.02222 

21201 = 

o.o ; 

¥(205= 

-0.33458 

P/PSI 20 5= 

0.02121 

-:«OO.NEWT 

.PRESS. AT 

INPUT COORD o 



r $m- = 

20 AMD XC 

5 5=' ■ 0* 

27419 



: ,TH€ Zb V 

AND P -COORDINATES ARE: 



• ’-21 .11= 

0.0 

¥8 1,1 = 

0.36120 

P/PSf 11 = 

0.34043 

-7*; 21= 

0.05945 

YC .25= 

0.35627 

P/PS? 21= 

0.33498 

•’ i21. 31= 

0.11728 

YC 3-1 = 

0.34163 

P/P-Sf 31 = 

0.32000 
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I 


I 


■I 


ZC 4)= 0.1=7191 
ZC 51- 0.22185 
ZC 61= 0.26574 
Z ( 71= 0.30238 
Zl 81= 0.33077 
ZC 91= 0.35014 
Z(10)= 0.35996 
ZC11)= 0.35996 
Z ( 12)= 0.35014 
ZC 13)= 0.33077 
Z ( 14)= 0.30238 
Z (15)= 0.26574 
ZC 16)= 0.22185 
Z ( 17)= 0.17191 
Z ( 18 1= 0.11728 
Z ( 19)= 0.05945 
Z C 20)= 0.0 

HOD. NEWT. PRESS. AT 


YC 4 )= 0.31766 

Yl 5)= 0.28503 

YC 61= 0.24463 

YC 7 )= 0.19756 

YC 8)= 0.14509 

YC 9) = 0.08867 

YC 10)= 0.02983 

YC 11 )= -0.02983 

YC 12 )= -0.08867 

YC 13)— -0.14509 

YC 14)= -0.19756 

YC 15)= -0.24463 

YC 16)= -0.28503 

YC 17 )= -0.31766 

Y ( 1 8 ) = -0.34163 

YC 19)= -0.35627 

YC 20)= -0.36120 

INPUT COORD. 


P/PSC 4)= 0.29659 
P/PSC 5)= 0.26581 
P/PSf 61= 0.22921 
P/PSC 7)= 0.19003 
P/PSC 8 ) = 0.15201 
P/PSC 91= 0.11558 
P/PSC10)= 0.08347 
P/PSC 11)= 0.05638 
P/PSC12 1— 0^03522 
P/PSC 13)= 0.02058 
P/PSC 14)= 0.01158 
P/PSC 15)= 0.00759 
P/PSC 16)= 0.00689 
P/PSC 17)= 0.00689 
P/PSC 18)= 0.00689 
P/PSC 19)= 0.00689 
P/PSC20)= 0.00689 


FOR JH = 20 AND XC 61= 0.28169 

THE Z»Y AND P COORDINATES ARE: 


ZC 1) = 

0.0 

YC 11= 

0.36321 

ZC 21= 

0.05978 

YC 2 ) = 

0.35825 

ZC 31= 

0.11793 

YC 31= 

0.34353 

ZC 41= 

0.17287 

VC 4 ) = 

0.31943 

ZC 5)= 

0-22309 

YC 5)= 

0.28662 

ZC 61= 

0.26722 

YC 61= 

0.24599 

ZC 7)= 

0.30406 

YC 71= 

0.19865 

ZC 8) = 

0.33261 

YC 8 ) = 

0.14590 

ZC 91= 

0.35209 

YC 9 ) = 

0.08916 

ZC10)= 

0.36197 

YC 101= 

0.02999 

ZCil)= 

0-36197 

YC 111= 

-0.02999 

ZC12)= 

0.35209 

YC 121= 

-0.08916 

Z C 1 3 ) = 

0.33261 

YC 13)= 

—0^14590 

ZC 141= 

0.30406 

YC 14)= 

-0.19865 

Z(15)= 

0.26722 

YC 151= 

—0^24599 


P/PSC 1)= 
P/PSC 21 = 
P/PSC 3)= 
P/PSC 41= 
P/PSC 5 ) = 
P/PSC 61 = 
P/PSC 7) = 
P/PSC 8 )= 
P/PSC 9) = 
P/PSC 101= 
P/PSC 111= 
P/PSC 121= 
P/PSC 13)= 
P/PSC 14)= 
P/PSC 15)= 


0.33039 

0.32499 

0.31020 

0.28718 

0.25685 

0.22071 

0.18219 

0.14503 

0.10943 

0.07830 

0.05215 

0.03198 

0.01833 

0.01027 

0.00715 
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21161* 

0.22309 

Y? 16 1= 

-0.28662 

P/PSI l6l= 

0.00689 

ZIL71= 

0.17287 

YC 17)= 

-0.31943 

P/PSC17I* 

0.00689 

ZI181= 

0. 11793 

YC 181= 

-0.34353 

P/PSf 131= 

0.00689 

Z«191= 

0.05978 

Yf 191= 

-0.35825 

P/PSI 191= 

0. 00689 

ZI201= 

0.0 

Yi 201= 

-0.36321 

P/PSI 201= 

0. 00689 

MOD. NEWT 

•PRESS. AT 

INPUT COORD. 




FOR JM - 

20 AND XI 

71= 

0.50000 



THE Z.Y 

AND P COORDINATES ARE s 



Zf li= 

0.0 

Yl 

1 1= 

0.42170 

P/PSf 11= 

0. 33439 

Zl 21 = 

0.06942 

YC 

21= 

0.41595 

P/PSf 21= 

0.32936 

u» 

II 

0.13693 

YC 

31 = 

0.39885 

P/PSI 31= 

0.31436 

Zl 41= 

0-20071 

YC 

41= 

0.37088 

P/PSl 41= 

0.29047 

Zl 51= 

0.25901 

VI 

51 = 

0.33278 

P/PSI 51= 

0.25960 

Zt 61= 

9.31026 

YC 

61 = 

0.28561 

P/PSI 61= 

0.22412 

Z( ,71 = 

0.35303 

YC 

71 = 

0.23065 

P/PSI 71= 

0.18612 

Zl 81= 

3.38618 

YC 

81 = 

0.16940 

P/PSI 81= 

0. 14780 

Zl 91 = 

0.40680 

YC 

91 = 

0.10352 

P/PSI 91= 

0.11181 

Zl 101= 

0.42026 

YC 101= 

0^03482 

P/PSI 101= 

0.07985 

Z 1111= 

0.42026 

Yl 111= 

-0.03482 

P/PSI111= 

0.05345 

ZI121= 

0.40880 

Y 1 121= 

-0.10352 

P/PSI 121= 

0.03324 

21131= 

0.38618 

YC 131= 

-0.16940 

P/PSI 131= 

0.01919 

21141= 

0.35303 

Yl 141= 

-0.23065 

P/PSI 141= 

0.0109 I 

21151= 

0.31026 

Yf 151= 

-0.28561 

P/PSC151= 

0.00730 

21161'= 

0.25901 , 

YC 161= 

-0.33278 

P/PSI 161= 

0* 00689 

21171= 

0.20071 

YC 171= 

-0.37088 

P/PSI 171= 

0.00689 

21181= 

0.13693 

YC 18 1= 

-0.39885 

P/PSI 181= 

0.00689 

Zl 191= 

0.06941 

Yl 191= 

-0.41595 

P/PSI 191= 

0.00689 

21201= 

0.0 

Yl 201= 

-0.42170 

P/PSI 201= 

0.00689 


NOD. NEWT. PRESS. AT INPUT COORD. 


FOR JR * 20 AND XI 8 1= 0.75000 

THE Z« V MO ? COORDINATES ARE s 
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z< n= 

0.0 

Yf n = 

0.48869 

P/PSf 1)= 

Zf 21= 

0.08044 

Yf 2) = 

0.48202 

P/PSf 21= 

Z< 31 = 

0.15868 

Yf 3) = 

0.46221 

P/PSf 3 ) = 

Z( 41 = 

0.23259 

Yf 41= 

0.42979 

P/PSf 41= 

Zf 51 = 

0.30016 

Yf 5 )= 

0.38564 

P/PSf 51= 

Zf .61= 

0.35954 

Yf 6 ) = 

0.33098 

P/PSf 61= 

Zf 71 = 

0.40911 

Yf 7) = 

0-26729 

P/PSf 71= 

Zf 8»= 

0.44753 

Yf 8 ) = 

0.19630 

P/PSf 81= 

Z( 9 ) = 

0.47373 

Yf 9) = 

0.11997 

P/PSf 91= 

ZflOI= 

0.48702 

Y f 1 0 1 = 

0.04036 

P/PSflOl= 

Zfll)= 

0.48702 

Yf 11)= 

-0.04036 

P/PSfll)= 

Z ( 12)= 

0.47373 

Yf 121= 

-0.11997 

P/PSf 121= 

Zf 131= 

0.44753 

Yf 13)= 

-0.19630 

P/PSf 131= 

Zf 141= 

0.40911 

Yf 14)= 

-0.26729 

P/PSf 14) = 

Z f 15 )= 

0.35954 

Yf 15)= 

-0.33098 

P/PSf 151= 

Z f 16) = 

0.30016 

Yf 16)= 

-0.38564 

P/PSf 161= 

Zf 17)= 

0.23259 

Yf 178= 

-0.42979 

P/PSf 17) = 

Z f 1 8 1 = 

0.15868 

Yf 188= 

-0.46221 

P/PSf 181= 

Z f 19) = 

0.08044 

Yf 191= 

-0.48202 

P/PSf 19)= 

Zf 205= 

0.0 

Yf 201= 

-0.48869 

P/PSf 201= 

MOD.NEWT 

.PRESS. AT 

INPUT COORD. 


FOR JM = 

20 AMO Xt 

91= 1. 

00000 


THE Z* Y 

AND P COORDINATES ARE; 


Zf 11 = 

0.0 

Yf 11= 

0.55568 

P/PSf 11= 

Zf 21 = 

0.09146 

Yf 25 = 

0.54810 

P/PSf 21= 

Zf 3S = 

0.18043 

Yf 3 ) = 

0.52557 

P/PSf 3 ) = 

Zf 4i = 

0.26447 

Yf 41= 

0.48870 

P/PSf 41= 

Zf 5 )= 

0.34130 

Yf 5 ) = 

0.43851 

P/PSf 5 ) = 

Zf 61 = 

0.40882 

Yf 61= 

0.37635 

P/PSf 61= 

Zf 7 )= 

0.46519 

Yf 71 = 

0.30393 

P/PSf 71= 

Zf 8 ) = 

0.50887 

Yf 8) = 

0.22321 

P/PSf 8 ) = 

Zf 91 = 

0.53867 

Yf 9 ) = 

0.13641 

P/PSf 9 ) = 

ZflO)= 

0.55378 

Yf 10)= 

0.04589 

P/PSflOl= 

Zf 11)= 

0.55378 

Yf 118= 

-0.04589 

P/PSflI)= 

Zf 121= 

0.53867 

Yf 12)= 

-0.13641 

P/PSf 12 1= 


0. 33344 
0.32830 
0.31336 
0.28964 
0.25895 
0. 22327 
0.18517 
0. 14709 
0.11122 
0.07950 
0.05313 
0.03293 
0.01898 
0.01074 
0.00726 
0.00689 
0.00689 
0.00689 
Ob 00689 
0.00689 


0.33369 
0.32859 
0.31361 
0.28988 
0. 25913 
0.22350 
0.18542 
0. 14726 
0.11140 
0.07959 
0.05321 
0.03301 
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Z I 13) = 

0.50887 

Y( 13)= 

-0.22321 

P/PSC 131= 

0,01903 

Zf 143 = 

0.46519 

YC 14) = 

-0.30393 

P/PSI 14) = 

0.01079 

ZI 151 = 

0.40882 

YC 153= 

-0.37635 

P/PSI 15) = 

0.00727 

ZIifc)= 

0.34130 

¥5 161= 

-0.43851 

P/PSI 163= 

0.00689 

m?i= 

G. 26447 

Y(17)= 

-0.48870 

P/PSI 171 = 

0.00689 

ZC 181= 

0.18043 

¥(183= 

-0.52557 

P/PSI 18)= 

0.00689 

ZI19!= 

0.0914& 

YC 191= 

-0.54810 

P/PSI 19 J= 

0.00689 

Z!20i= 

0.0 

Y( 20 1= 

-0.55568 

P/PSC 203= 

C. 00689 

HOD. NEWT 

.PRESS. AT 

INPUT COORD. 




FOR JM = 

‘ 20 AND XI 

103= 2. 

00000 



THE Z«Y 

AND P COORDINATES ARE: 



Z( 11= 

0.0 

¥f 1> = 

0.82363 

P/PSC 11= 

0.33360 

Zf 23= 

0.13556 

Y { 2 ? = 

0.81239 

P/PSC 23= 

0.32848 

ZI 33 = 

0.26743 

YC 33 = 

0.77900 

P/PSC 31= 

0.31353 

It 4i = 

©.'39200 

Y( 4 3 = 

0.72436 

P/PSI 41= 

0.28982 

ZI 5)= 

0.50588 

Y( 5 ) = 

0.64996 

P/PSI 53= 

0.25907 

ZI 63= 

£3.60596 

YC 63 = 

0.55783 

P/PSI 61= 

■ 0.22343 

ZI 71= 

0.68951 

YC 7 S = 

0.45048 

P/PSI 73= 

0.18533 

ZI 81 = 

0.75425 

YC 83 = 

0.33085 

P/PSI 83= 

0.14722 

ZI 91 = 

0.79842 

YC 9 3 = 

0.20219 

P/PSI 93= 

0.11134 

21101= 

0.82081 

YC 103= 

0.06801 

P/PSUG) = 

0.07955 

21113= 

0.82081 

YC 113= 

-0.06801 

P/PSI 11 3= 

0.05319 

Z 1 12) = 

0.79842 

YC 123= 

-0.20219 

P/PSI 121= 

0.03299 

Z 1 133= 

0.75425 

Y C 1 3 3 = 

-0.33085 

P/PSI133= 

0.01902 

ZI 143= 

0.68951 

YC 143= 

-0.45048 

P/PSI 141= 

0.0107? 

21153= 

G. 60596 

YC 153= 

-0.55783 

P/PSI 151= 

0.00727 

Z 1 1 6 1= 

0,50588 

Y( 161= 

-0.64996 

P/PSI 163= 

0.00689 

Z 1173= 

0.39200 

Y! 173= 

-0.72436 

P/PSC173 = 

0.00689 

21181= 

0.26743 

YC 183= 

-0.77900 

P/PSI 181= 

0. 00689 

ZI19i= 

0.13556 

YC 193= 

-0.81239 

P/PSI 193= 

0.00689 

Z 1201= 

0.0 

Y C 203 = 

-0.82363 

P/PSI20i= 

0.00689 


HOD .NEWT .PRESS. AT INPUT COORD. 
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I 

I 


FOR JM = 20 AND XI 111= 3.00000 

THE 2» Y AND P COORDINATES AREs 


ZMJ* 0.0 
Zl 21= 0-17967 

Zl 31= 0.35443 

Zl 41 = 0.51953 

ZC 51= 0.67046 

Zl 61= 0.80310 

Zl 71= 0.91383 

Zl 81= 0.99963 

Zl 91 = 1.05817 

Z 1 101= 1.08785 

Zllll= 1.08785 
Z 1 121=1 1.05817 

Z C 1 31= 0.99963 

Z 1 141= 0.91383 

Z I 151= 0.80310 

Z I 16 1= 0.67046 

Z I 17 1 = 0.51953 

Z 1 181= 0.35443 

ZI191= 0-17967 

ZI201= 0.0 

MOD. NEWT .PRESS. AT 


Yl 11 = 

1.09157 

Yl 21 = 

1.07669 

Yl 31 = 

1.03243 

Yl 41 = 

0.96001 

Yl 51 = 

0.86141 

Yl 61= 

0.73930 

Yl 71 = 

0.59703 

Yl 81 = 

0.43848 

Yl 91 = 

0.26797 

Yl 101= 

0.09014 

Yl 11 1 = 

-0.09014 

Yl 121= 

-0.26797 

Yf 131= 

-0.43848 

Yl 141= 

-0.59703 

Yl 151= 

-0.73930 

Yl 161= 

-0.86141 

Yl 171= 

-0.96001 

Yl 181= 

-1.03243 

Yl 191= 

-1.07669 

Y 1201= 

-1.09157 


COORD. 


P/PS I 11= 
P/PSI 21= 
P/PS1 31= 
P/PSI 41= 
P/PSI 51= 
P/PSI 61= 
P/PSI 71= 
P/PSI 81= 
P/PSI 91= 
P/PSI 101= 
P/PSI111= 
P/PSI 12 1 = 
P/PSI 131= 
P/PSI 141= 
P/PSI 151= 
P/PSI 161= 
P/PSI 171= 
P/PSI 181= 
P/PSI 191= 
P/PSI20I= 


0.33361 
0.32851 
0.31354 
0.28983 
0.25908 
0.22343 
0.18535 
0.14723 
0.11135 
0.07956 
0.05319 
0. 03299 
0.01902 
0.01077 
0.00727 
0.00689 
0.00689 
0.00689 
0.00689 
0.00689 


FOR JM = 20 AND XI 121= 4.00000 

THE Z„Y AND P COORDINATES AREs 


ZC 11= 0.0 YC 11= 1.35952 P/PSI 11= 0.33362 
Zl 21= 0.22377 Yl 21= 1.34098 P/PSI 21= 0.32849 
Zl ,31= 0.44144 Yf 31= 1.28586 P/PSI 31= 0.31354 
Zl 41= 0.64706 Yl 41= 1.19567 P/PSI 41= 0.28983 
Zl 51= 0.83504 Yl 51= 1.07286 P/PSI 51= 0.25907 
ZC 61= 1.00023 YC 61= 0.92078 P/PSI 61= 0.22344 
Zl 71= 1.13815 Yl 71= 0.74359 P/PSI 71= 0.18535 
ZC 81= 1.24502 Yl 81= 0.54611 P/PSI 81= 0.14723 
Zl 91= 1.31792 Yl 91= 0.33374 P/PSI 91= 0.11134 
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ZU0J= 

1,35488 

YC 101= 

0.11227 

P/PSI 101= 

0.07955 

Zflll- 

1-35488 

YC 111= 

-0.11227 

P/PSC 111 = 

0-05319 

Z ( 121= 

1-31792 

YC 121= 

-0.33374 

P/PSC 121 = 

0.03299 

Z 1 131= 

1-24502 

YC 131= 

-0.54611 

P/PS! 131= 

0.01902 

Z C 141= 

1,13815 

Yf 141= 

-0.74359 

P/PSC 141= 

0.01077 

Zf 151= 

1-00023 

Y 1-151= 

-0.92078 

P/PSC 151= 

0.00727 

Z C 161= 

0.83504 

YC 161= 

-1.07286 

P/PSC 161= 

0.00689 

Z C17»= 

0-64706 

YC 17 1 = 

-1.19567 

P/PSI 171= 

0.00689 

ZC18S= 

0-44144 

YC 181= 

-1.28586 

P/PSC 181= 

0.00689 

Z C 191 = 

0.22377 

YC 191= 

-1.34098 

P/PSC 191 = 

0.00689 

Z C 201= 

0.0 

YC 201 = 

-1.35952 

P/PSC 201= 

0.00689 

MOD. NEWT 

.PRESS. AT 

INPUT COORD- 




FOR JM 

= 20 AND XC 131= 5. 

50000 



THE Z, i 

? AND P COORDINATES ARE: 



ZC 11= 

0-0 

YC 11= 

1.76145 

P/PSC 11= 

0.33361 

Zf 21= 

0.28992 

YC 21 = 

1.T73742 

P/PSC 21= 

0.32850 

ZC 31= 

0.57194 

YC 31 = 

1.66601 

P/PSC 31= 

0.31354 

ZC 41 = 

0.83836 

YC 41 = 

1.54915 

P/PSC 41= 

0.28982 

Zf 51 = 

1.08190 

YC 51 = 

1.39003 

P/PSC 51= 

0.25907 

Za 61 = 

1.29594 

YC 61 = 

1.19300 

P/PSC 61= 

0.22344 

ZC 71 = 

1.47462 

YC 71 = 

0-96342 

P/PSC 71= 

0.18535 

ZC 81= 

i. 61309 

YC 81 = 

0.70757 

P/PSC 81= 

0.14724 

ZC 91 = 

1.70755 

YC 91 = 

0.43241 

P/PSC 91= 

0.11134 

ZC 101= 

1.75543 

YC 101= 

0.14546 

P/PSC101= 

0.07955 

ZC111= 

1.75543 

YC 111= 

-0.14546 

P/PSC11)= 

0.05320 

Z C 1 2 1 = 

1.70755 

YC 121= 

-0.43241 

P/PSC 121= 

0.03299 

Z C 131 = 

1.61309 

Y C 1 3 1 = 

-0.70757 

P/PSC 131= 

0.01902 

Z C 141= 

1.47462 

YC 141= 

-0.96342 

P/PSC 141= 

0.01077 

Z 5 1 5 1 = 

1.29594 

YC 15 1= 

-1.19300 

P/PSI 151= 

0.00727 

Z 1 161= 

1.08190 

YC 16 1= 

-1.39003 

P/PSC 161= 

0.00689 

Z C 17 } = 

0.83836 

YC 17 1= 

-1.54915 

P/PS f 171= 

0.00689 

Z C 18 1 = 

0.57194 

YC 18 1= 

-1.66601 

P/PSC 181 = 

0.00689 

Z C 191 = 

0.28992 

YC 191= 

-1.73742 

P/PSC19) = 

0.00689 

Z«201= 

0.0 

YC 20 1= 

-1.76145 

P/PSC201= 

0.00689 
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MOD.NEWT. PRESS. AT INPUT COORO. 

FOR JM = 20 AND XI 141= 7.00000 

THE Z.Y ANO P COORDINATES ARE: 


ZC 11 = 

0.0 

ZC 21 = 

0.35608 

ZC 31 = 

0.70245 

ZC 41 = 

1.02965 

ZC 51 = 

1.32877 

ZC 61 = 

1.59164 

ZC 71 = 

1.81110 

ZC 81= 

1.98116 

ZC 91 = 

2.09717 

ZC101= 

2.15598 

ZC111= 

2.15598 

Z C 121= 

2.09717 

ZC131= 

1.98116 

ZC 141= 

1.81110 

Z C 15 1= 

1.59164 

Z C 161= 

1.32877 

Z C 171= 

1.02965 

ZC 181= 

0.70245 

Z ( 191= 

0.35608 

Z (201= 

0.0 

MOO. NEWT 

-PRESS. AT 


YC 11= 

2.16337 

Y( 21 = 

2.13387 

YC 31= 

2.04615 

YC 41 = 

1.90263 

YC 51 = 

1.70720 

YC 61 = 

1.46521 

YC 71 = 

1.18325 

YC 81 = 

0.86902 

YC 91 = 

0.53108 

YC i0i= 

0.17865 

YC 111= 

-0.17865 

YC 121= 

-0.53108 

YC 131= 

-0.86902 

YC 141= 

-1.18325 

YC 151= 

-1.46521 

YC 161= 

-1.70720 

YC 171= 

-1.90263 

YC 181= 

-2.04615 

YC 191= 

-2.13387 

Y C 201 = 

-2.16337 


INPUT COORD* 


P/PSC 11 = 0.33361 
P/PSC 21= 0.32850 
P/PSC 31= 0.31354 
P/PSC 4)= 0.28982 
P/PSC 51 = 0.25907 
P/PSC 61= 0.22344 
P/PSC 71= 0.18535 
P/PSC 81= 0.14723 
P/PSC 91= 0.11134 
P/PSCiOl- 0.07955 
P/PSC 111= 0.05320 
P/PSC 121= 0.03299 
P/PSC 131= 0.01902 
P/PSC 141= 0.01077 
P/PSC 151 = 0.00727 
P/PSC 161= 0.00689 
P/PSC 171= 0.00689 
P/PSC 181= 0.00689 
P/PSC 191= 0.00689 
P/PSC 201= 0.00689 


FOR JM = 20 ANO XC151= 10.00000 

THE Z.Y AND P COORDINATES ARE: 


ZC 

11= 

0.0 

YC 

11 = 

ZC 

21 = 

0.48839 

YC 

21 = 

ZC 

31 = 

0.96345 

YC 

31= 

ZC 

41 = 

1.41224 

YC 

41 = 

ZC 

51= 

1.82250 

YC 

51= 

ZC 

61= 

2.18305 

YC 

61= 


2.96722 

P/PSC 11= 

0.33362 

2.92675 

P/PSC 21= 

0.32850 

2.80645 

P/PSC 31= 

0.31354 

2.60959 

P/PSC 41= 

0.28982 

2.34155 

P/PSC 51=; 

0.25907 

2.00964 

P/PSC 61= 

0.22344 
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ZC 7) = 

2.48406 

YC 7 ) = 

1.62291 

P/PSC 7)= 

0.18535 

ZC 81= 

2.71730 

YC 81 = 

1.19192 

P/PSC 81= 

0.14723 

Zl 91 = 

2.87642 

YC 91 = 

0.72841 

P/PSC 91= 

0.11134 

zcioi=: 

2.95708 

YC 101= 

0.24503 

P/PSC10J= 

0.07955 

ZC11)= 

2-95708 

Yi in-: 

-0.24503 

P/PSClll- 

0.05320 

Z 1 12)= 

2.87642 

YC 121= 

-0.72841 

P/P SCI 21= 

0.03299 

£1131= 

2.71730 

YC 131= 

-1.19192 

P/PSC 131= 

0.01902 

Z C 1 4 S = 

2.48406 

YC 14 1= 

-1.62291 

P/PSC 141= 

0.01077 

21151= 

2.18305 

YC 151= 

-2.00964 

P/PSS 15B= 

0.00727 

ZS16) = 

i. 82250 

YI 16)= 

-2.34155 

P/PSC 16?= 

0.00689 

ZI17?= 

1.41224 

YC 17)= 

-2.60959 

P/PSC 173= 

0.00689 

mss= 

0.96345 

¥C 181= 

-2.80645 

P/P SI 18)= 

0.006S9 

1119)= 

0.48839 

YC 191= 

-2.92675 

P/PSC 19)= 

G. 00689 

Z i 201= 

0.0 

Yf 20)= 

-2.96722 

P/PSC2G3= 

0.00689 

NGD.NEiir 

.PRESS. .41 

INPUT COORD. 




FOR JM = 

20 AND M 

161= 13. 

00000 



THE Z.Y 

AND. P COOR 

DINATES ARE; 



Zi 11= 

0.0 

YC 11 = 

3.77107 

P/PSC 11= 

0.33361 

Zl 21 = 

0.62070 

YC 21 = 

3.71963 

P/PSC 21= 

0.32850 

ZI 31= 

1.22446 

Y( 31 = 

3.56674 

P/PSC 31= 

0.31354 

Zl 4) = 

1.79483 

YC 4 ) = 

3.31655 

P/PSC 41= 

0.28983 

ZC 5 ) = 

2.31624 

YC 51 = 

2.97590 

P/PSC 51= 

0.25907 

Zl 6 3 = 

2.77446 

Yi 6 3= 

2.55407 

P/PSC 63= 

0.22344 

Zl -71= 

3.15701 

YC 71 = 

2.06258 

P/PSC 71= 

0.18535 

Z! 83 = 

3.45344 

YC 81 = 

1.51482 

P/PSC 3 ) = 

0.14723 

ZC 9 1 =: 

3.65567 

YC 9) = 

0.92574 

P/PSC 91= 

0. 11134 

Zf 10)= 

3.75819 

YC 101= 

0.31141 

P/PSC 101= 

0.07955 

ZC113= 

3.75819 

YC 111— 

-0.31141 

P/PSC1I1= 

0.05320 

ZC12)= 

3.65567 

YC 12)= 

-0.92574 

P/PSC 121= 

0.03299 

ZC13)= 

3.45344 

YC 13)= 

-1.51482 

P/PSC 131= 

0.01902 

Zf 143= 

3.15701 

YC 14)= 

-2.06258 

P/PSC 14) = 

0.01077 

Z C 151= 

2.77446 

Y(15>= 

-2.55407 

P/PSC 153= 

0.00727 

Z S 161 = 

2.31624 

YI 16)= 

-2.97590 

P/PSC 161= 

0.00689 

ZC 171= 

1.79483 

Y(17)= 

-3.31655 

P/PSC 171= 

0.00689 

Z? 181= 

1.22446 

YC 18)= 

-3.56674 

P/PSC 18 ) = 

0.00689 
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Z ( 19 1 = 0.62070 YC 19 1= -3.719j63 

Z f 201= 0.0 YC 20)= -3.77107 

MOD. NEWT. PRESS. AT INPUT COORD. 


P/PSC19) 
P/PSC 201 : 


FOR JM = 

20 ANO XC 17 1= 16. 

00000 


THE Z* Y 

AND > P COORDINATES ARE: 


Zt 11 = 

0.0 

YC 11= 

4.57491 

P/PSC 11= 

ZC 21= 

0.75301 

YC 21 = 

4.51252 

P/PSC 21= 

ZC 31= 

1.48547 

YC 31= 

4.32703 

P/PSC 31= 

ZC 41 = 

2.17742 

YC 41 = 

4.02352 

P/PSC 41= 

ZC 51 = 

2.80997 

YC 51= 

3.61025 

P/PSC 51= 

ZC 61= 

3.36587 

YC 61 = 

3.09850 

P/PSC 6 1 — 

ZC 71 = 

3.82996 

YC 71 = 

2.50224 

P/PSC 71= 

Zt 81= 

4.18958 

YC 81 = 

1.83772 

P/PSC 81 = 

ZC 91 = 

4.43492 

YC 91 = 

1.12307 

P/PSC 91= 

ZC101= 

4.55929 

YC 10 ) = 

0.37779 

P/PSC 101= 

ZClli= 

4.55929 

YC 111= 

-0.37779 

P/PSC1U = 

Z C 1 2 1 = 

4.43492 

YE 121= 

-1-12307 

P/PSC 121= 

Z 1 131 = 

4.18958 

YC 131= 

-1.83772 

P/PSC 131= 

Zt 141= 

3.82996 

YC 141= 

-2.50224 

P/PSC 141= 

ZC 151= 

3.36587 

YC 151= 

-3.09850 

P/PSC151 = 

Z C 161= 

2.80997 

YC 161= 

-3.61025 

P/PSC 161= 

Z C 171= 

2. 17742 

YC 171= 

-4.02352 

P/PSC 171= 

Z 1 18 1 = 

1-48547 

VC 181 = 

-4.32703 

P/PSC 181= 

Z C 191 = 

0.75301 

YC 191= 

-4.51252 

P/PSC 191= 

Z C 20 1 = 

0.0 

YC 201 = 

-4.57491 

P/PS I 201= 

MOO. NEWT 

-PRESS. AT 

INPUT COORD. 


FOR JM = 

20 AND XC 181= 19. 

00000 


THE Z,Y 

AND P COORDINATES ARE: 


ZC 11= 

0.0 

YC 11 = 

5.37876 

P/PSC 11= 

ZC 21= 

0.88532 

YC 21 = 

5.30540 

P/PSC 21= 

ZC 31= 

1.74648 

YC 31= 

5.08733 

P/PSC 31= 


0.00689 

0.00689 


0.33361 
0.32650 
0.31354 
0.28983 
0.25907 
0.22344 
0.18535 
0. 14723 
0.11134 
0.07955 
0.05319 
O'. 03299 
0.01902 
0.01077 
0.00727 
0.00689 
0.00689 
0.00689 
0.00689 
0.00689 


0.33362 

0.32850 

0.31354 


213 


N N M N N N (M H N N N 


ZI 4) = 

2.56001 

Y( 4) = 

4.73048 

P/PSi 41= 

ZI 51= 

3.30370 

Y( 5 1= 

4.24460 

P/P SI 51= 

Zl 61= 

3.95728 

Y( 6 ) = 

3.64294 

P/PSC 61= 

ZI 7) = 

4.50292 

Yl 71 = 

2.94190 

P/PSI 7) = 

Zi 81 = 

4.92573 

Yl 81 = 

2.16062 

P/PSI 8)= 

ZC 9)= 

5.21417 

YI 9} = 

1.32041 

P/PSI .'91 = 

ZU0J= 

5.36039 

YC 10J= 

0.44417 

P/PS(10}= 

ZC1I)= 

5.36039 

Y( 11) = 

-0.44417 

P/PSI IS ) = 

Z 1121= 

5.2141? 

Y( 12 1 = 

-1.32041 

P/PSI 12) = 

ZI 131= 

4.92573 

Y( 13) = 

-2.16062 

P/PSC 13)= 

ZI 141= 

4.50292 

Y( 141 = 

-2.94190 

P/PSI 141= 

Z'{ 151* 

3.95728 

YI 15) = 

-3.64294 

P/PSI 155 = 

Z I 161= 

3-30370 

Yt 161= 

-4.24460 

P/PSI 161= 

It 17}= 

2.56001 

YI 17) = 

-4.73048 

P/PSI 171= 

21183= 

1.74648 

¥118)= 

-5.08733 

P/PSI 18>= 

Zt 19)= 

0.88532 

YC 19)= 

-5.30540 

P/PSI 191= 

Z 1 20}= 

0.0 

YC20)= 

-5.37876 

P/PSI 201= 


MOO. NEWT ® PRESS . AT INPUT COORD-i 


FOR JM - 20 AND Mfi9l= 22.00000 
THE Z*¥ &ND P COORDINATES ARE: 


ZI U= 

0.0 

YI 

15= 

6.18261 

P/PSC 11= 

ZI 2) = 

1.01762 

YC 

23= 

6.09829 

P/PSC 21= 

ZI 3)= 

2.00749 

YI 

3 ) = 

5.84762 

P/PSC 33= 

ZI 41 = 

2-94260 

Yl 

4 3 = 

5.43744 

P/PSC 4) = 

ZI 55 = 

3.79744 

Y? 

53= 

4.87895 

P/PSC 5 3 = 

ZI 63 = 

4.54869 

Y( 

6 1= 

4.18737 

P/PSC 61 = 

ZI 7) = 

5.17587 

YC 

71 = 

3.38157 

P/PSC 71= 

ZI 85 = 

5-86187 

YC 

83 = 

2.48353 

P/PSI 8 )= 

ZI 95 = 

5.99342 

YC 

9> = 

1.51774 

P/PSI 9)= 

ZI 101= 

6.1,6149 

Yl 10)= 

0.51056 

P/PSI10I= 

2111)= 

6.18149 

YC 11)= 

-0.51056 

P/PSfll)= 

ZC12)= 

5.99342 

Y(12)= 

-1.51774 

P/PSI 121= 

Z I 13)= 

5.66187 

Y C 1 3 1 = 

-2.48353 

P/PSC 131= 

*** 

>— 

* 

II 

5-17587 

YC 14)= 

-3.38157 

P/PSI 14)= 

Ztl5)= 

4.54869 

Y(i5)= 

-4.18737 

P/P SI 15 ) = 


0.28982 
0.25907 
0. 22344 
0.18535 
0.14723 
0.11134 
0.07955 
0.05320 
0.03299 
0.01902 
0.01077 
0.00727 
0.00689 
0.00689 
0.09689 
0.00689 
0.00689 


0.33361 

0.32850 

0.31354 

0.28983 

0.25907 

0.22344 

0.18535 

0.14723 

0.11134 

0.07955 

0.05320 

0.03299 

0.01902 

0.01077 

0.00727 
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Z ( 161= 
Z < 17 1= 
Z(181= 
Z ( 19 1= 
Z (201= 

3.79744 

2.94260 

2.00749 

1.01762 

0.0 

Y( 161= 
Y( 17 1= 
Y( 18 1 = 
Y( 191= 
Y( 201= 

-4.87895 

-5.43744 

-5.84762 

-6.09829 

-6.18261 

P/PS (161= 
P/PSC 17 1= 
P/PSf 181= 
P/PSC19)= 
P/PS(20)= 

0.00689 
0.00689 
0.00689 
0.00689 
0. 00689 

STREAMLINE 

CALCULATIONS 





INPUT DATA: 

ANG. OF ATTACK 3 0.20000E 02 DEG. 
ZETA(WALL)= 0.25100E 00 

NO. OF STREAMLINES TO BE CALCULATED = 2 


B = Q.99335E 00 

STAGNATION PQINT»XO = 0.22615E-01 

TW = 0. 530067E 03 DEG. R QWS = 0.359411E 02 BTU/SO FT-SE 

HSP = 0.227223E-01 8TU/SQ FT-SEC-DEG R 


FOR STREAMLINE NO. 1 WHERE BETA = 0.10000E 01 DEG. 
THE CALCULATIONS ARE AS FOLLOWS: 


X 

Y 

Z 

H 

S 

BBAR 

P/PS 

QW/QWS 

UE 

HE 

RHOE 

HUE 

MACH NO 

RN/FT 

MOM THK 

RNM 

RNl 

R 

PHI 

DR/DPHI 

D2R/DPHI2 

OR/DX 

D2R/DX2 

QB/QSB 

0.37092E— 01 0. 

16261E 00 0 

-65 OOZE— 03 0 

• 37250E— 01 

0. 37256E— 01 0 

•68442E 00 

0.99188E 00 

0.91683E 00 

0.24302E 03 

0.12654E 08 

0. 10593E— 03 

0.95220E— 06 

0.10802E 00 

0.27036E 05 

0.14879E— 03 

0.40228E 01 

0.93173E— 02 

0.16261E 00 

0.22905E 00 

-0. 35879E— 05 

-0.9611 3E-04 

0.20780E 01 

— 0. 32710E 02 

0.91726E 00 
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0.38723E-01 0.I6596E 00 0. 

0.4096 IE— 01 Q» 62739E 00 0 


0.27279E 03 Q.12646E 08 

Q.L2129E 00 0.30313E 05 

0.8912IE 01 G.16596E 00 

—0.939406—04 0.2G262E 01 


71488E-03 0 

.98977E 00 
0. 10577E— 03 
0.144076-03 
0.24681E 00 
— 0.30774E 02 


0.42C86E— G1 C.&7260E GO 0.84333E-03 0 

0.484316-01 0.57318E 00 0.9B498E 00 

0-3308 IE 03 0.12629E 08 0.10541E-03 

0. 1471.96 00 0. 3666 3E- 05 0.13969E-03 

0.29715E 02 0.17260E 00 0.27994E 00 

— 0.58T45E— 04 0.19287E 01 -0.27364E 02 


0.49205E-01 0.IB569E 00 0.IO972E-02 0. 

0.83331 E— 01 0.54 11 IE 00 0.97310E 00 

0 « 44369 E 03 0.12585E 08 0.10450E-03 

0.197756 OS 0.488486 05 0.13759E-03 

0.828616 02 0.18569E 00 0.338546 00 

— 0.74779E— 04 G.17545E 01 -0.21995E 02 

0.5 683 8 £—0 1 0.I9848E 00 0.I3481E-02 0 

0.78233.6-01 0.53646E 00 0.95817E 00 

0.5 5460 E 03 0.12530E 08 0.1Q335E-03 

0.24782E 00 D.60568E 05 0.13774E-03 

0.250836 03 0.19849E 00 0.38916E 00 

— 0.S5729E— 04 0.16029E 01 -0.18031E 02 


0.649736-01 0.21096E 00 0.15964E-02 0 

0.931 3 IE— Cl Oc 54053E 00 0.94023E 00 

0. 66544 £ 03 0.12462E 08 0.101966-03 

0.298056 00 0.71S99E 05 0. 138516-03 

0.253 HE 03 0.21097E 00 0.43355E 00 

—G. 312106— 04 0.146966 01 -0.150396 02 


0. £26986-01 
0.122936 00 
0.886686 03 


0.234906 00 
0.554056 00 
0. 1 22916 08 


0.208386-02 
0.89566E 00 
0.984886-04 


40984E-G1 
S.94247E 00 
0. 95 186 £—06 
0. 43671 E 01 
—Q. 378986— 05 
0.94303E 0© 

» 483776-01 
5.966546 00 
0.951096-06 
0.51216E 01 
-0.42541 E-05 
2.967376 00 

, 629746—01 
3.974926 00 
0.94915E— 06 
0.6721 IE 01 
-0.486206-05 
0.976436 00 

.773926-01 
3.968546 00 
0.94670E— 06 
0.834256 01 
— Q. 52737E-05 
0.97090E 03 

» 91650E— 01 
G.95743E 00 
0. 943706-06 
0.99588E 01 
— 0.51271E— 05 
0.96079E OG 

. 11964E 00 
0.92620E 00 
0.93603E— 06 
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0.39990E 00 
0. 43663 E 03 
0.33884E-04 


0.93295E 05 
Q.23491E 00 
0.12443E 01 


0.14105E-03 
0.50826E 00 
-0.1092 IE 02 


0.13160E 02 
— 0. 35900E— 05 
0.93198E 00 


0.1022TE 00 0.25736E 00 0.25574E-02 0.14684E 00 

0.15273E 00 0.56701E 00 0.84157E 00 0.88306E 00 

0.1I045E 04 0.12G74E 08 0.94202E-04 0.92627E-06 

0.50258E 00 0. 11232E 06 0.14508E-03 0.16296E 02 

0.69451E 03 0.25737E 00 0.56933E 00 -0.40027E-06 

0.10582E-03 0.10596E 01 -^0.82818E 01 0.89165E 00 


0.14644E 00 0.29727E 00 0.34555E-02 0.1 984 IE 00 

0.21233E 00 0.61310E 00 0-71573E 00: 0.77453E 00 

0.15205E 04 0.11528E 08 0.83910E— 04 0.90I25E-06 

0.70808E 00 0.14156E 06 0. 15741 E-03 0.22283E 02 

0.1330IE 04 0.29729E 00 0.66599E 00 0.65593E-05 

0.15928E— 03 0-76852E 00 — 0.51671E 01 G.78893E 00 

0.19637E 00 0.32969E 00 0.42718E-02 0.24528E00 

0.27193E 00 0.85924E 00 0.56690E 00 0.6562IE 00 

0.19487E 04 G.10785E 08 0.71039E-04 0.86613E-06 

0»93823E 00 0.15983E 06 0.17505E— 03 0.27978E 02 

0.20827E 04 0.32972E 00 0.74234E 00 0.70036E-05 

0.11879E— 03 0.54246E 00 — 0.4I477E 01 0.67649E 00 

0.20960E 00 0.33653E 00 0.44595E-02 0.25606E 00 

0.28683E 00 0.97469E 00 0.52278E 00 0.62243E 00 

0.20715E 04 0.10538E 08 0.67045E-04 0.85417E-06 

0.10089E 01. 0.16259E 06 0.18096E-03 0.29422E 02 

0.22830E 04 0.33656E 00 0.75 92 0E 00 0. 76488 E— 05 

0.15989E— 03 0.49264E 00 -0.41644E 01 0.64425E 00 

0.23683E 00 0.34862E 00 0.48U6E-02 0.27627E 00 

0.316636 00 0.10888E 01 0.43U9E 00 0.54209E 00 

0.23280E 04 0.99738E 07 0.58427E-04 0.82624E-06 

0.11655E 01 0.16462E 06 0.19650E-03 0.32349E02 

0.26903E 04 0.34866E 00 0.79074E 00 0.41804E-05 

0.27427E— 03 0.39713E 00 -0.37105E 01 0.56632E 00 


217 


0.26494E 00 
0.34643E 00 
0.2 5431 E 04 
0.1 308 4£ 01 


0.35848E 00 
0.88841E 00 
0.94485E 07 
0. 16238E 06 


0.51324E— 02 
0.35680E 00 
0.51035E-04 
0.21602E— 03 


0 . 


0.309696 04 
0.13140E— 03 


0.35851E 00 
Q.30505E 00 


0.820276 00 
—0.237276 01 


O.27209E 00 0.36058E 00 0.52080E-02 ©« 

0.35388E 00 ' C.15869E 0& 0.34369E OO C 

0.25829E 04 0-93480E 07 O.49689E-04 

0.13357E 01 0. 16 1596 06 C.22G77E-03 

0.31975E 04 Q.36061E 00 0.82750E 00 

O. L9104E— 04 0.28238E 00 -0.18957601 

0.273.89E 00 0.381086 00 0.52267E-02 0. 

0.35574E 00 G.72013E 00 0.34085E 00 i 

0.25914E 04 0.93258E 07 0.49395E-04 

0.13417E 01 0. 161405 06 0.22191E-03 

0.32226E 04 0.36112E 00 G.82931E 00 

— 0. 17431 E— 04 0.27672E OO -0-17674E 01 

0.27748E 00 0.38205E 00 0.52637E-02 0. 

0.35947E 00 0.64877E 00 0.33569E 00 < 

0.26070E 04 0.92853E 07 0.48860E-04 

0.13528E01 0.L6105E 06 0.22410E-03 

0.32726E 04 0.36209E 00 0.83293E 00 

-0.I93398E-04 0.26775E 00 -0.15420E 01 

0.28 LOSE 1 00 0.36301E GO 0.53005E-02 0, 

0.36319E 00 ©.600156 00 0.33109E 00 I 

0.26210E 04 0.92488E 07 0.48380E-04 

G.13627E 01 0.16071E 06 0.22619E-03 

G.33226E 04 0.36305E 00 0.83655E 00 

— 0. 18377E— 03 Q.26428E 00 -0.14015E 01 

G. 282 886 00 0.36349E00 0.53189E-02 0 

0.365066 00 0.583506 00 0.32893E 00 

0.26276E 04 0.92315E 07 0.48155E-04 

0-13674E 01 0.16O55E 06 0.22721E-03 


29469E 00 
&.45950E 00 
0.79943E-06 
0-35077E 02 
— 0.618 19E— 05 
0.48423E 00 

>299026 00 
3.44087E 00 
0.79422E-06 
0.356756 02 
-0.10934E-04 
0.46538E 00 

. 30009E 00 
3.43638E 00 
0. 79306 E-06 
0.35817E 02 
-0. 121 30E—04 
Q.46081E 00 

» 30221 E 00 
C.42796E 00 
0. 79094E-06 
0.36090E 02 
-0.14613E-04 
0.45223E 00 

. 30433E 00 
0.42068E 00 
0.78903 E-06 
0.36352E 02 
—0. 17388E— 04 
0.44480E 00 

.305386 00 
0.41741E 00 
0.78813E— 06 
0.36478E 02 


218 


0.33475E 04 0.36352E 00 0.83836E 00 -0.18345E-04 

-0.22556E-03 0.26432E 00 -0.13586E 01 0.44148E 00 

0.2S648E 00 G.36444E 00 0.53559E-02 0. 3075 IE 00: 

0.36878E 00 0.55080E 00 0.32487E 00 0.41115E 00 

0.26400E 04 0.91988E 07 0.47729E-04 0. 78641 E-06 

0.13763E 01 „ 0.16023E 06 0.22919E-03 0.36724E 02 

0-33973E 04 0.36448E 00 0.84199E 00 — 0.22485E-04 

-0.30523E-03 0.26454E 00 -0.12767E 01 0.43509E 00 

0.29368E 00 0.36634E 00 0.54305E-02 0.31178E 00 

0. 37623E 00 0.48516E 00 0.31765E 00 0.39954E 00 

0.26622E 04 0.91399E 07 0.46969E-04 0.78331E-06 

0.13923E 01 0.15963E 06 0.23294E-03 0.37186E 02 

G.34966E 04 0.36638E 00 0.84927E 00 -0.26281E-04 

— 0.44934E— 03 0.26497E 00 -0. 11191E 01 0.42322E 00 

0.30809E 00 0.37016E 00 0.55820E-02 0.32047E 00 

G.39113E 00 0.35520E 00 0.30651E 00 0.37961E 00 

0.26968E 04 0.90472E 07 0.45787E-04 0.77840E-06 

0. 14176E 01 0.1S863E 06 0.23955E-03 0.38000E 02 

Q.36941E 04 0.37020E 00 0.86395E 00 -0.33537E-04 

-0.68059E— 03 0.26575E 00 -0.82829E 00 0.40273E 00 

0.33688E 00 0.37783E 00 0.58956E-02 0.33844E 00 

0.42093E 00 0.11S40E 00 0.29500E 00 0.34963E 00 

0.27331E 04 0.89488E 07 0.44552E-04 0.77317E-06 

0.14446E 01 0. 15749E 06 0.24922E^03 0-39248E 02 

0.40864E 04 0.3778SE 00 0.89396E 00 -0.410I1E— 04 

-0.94537E— 03 0«>26705E 00 -0.34105E 00 0.37154E 00 

0.39445E 00 0.39325E 00 0.65726E-02 0.37724E 00 

0.48053E 00 0.0 0.29956E 00 0.32835E 00 

0.27186E 04 0.89881E07 0.45043E-04 0.77526E-06 

0.14338E 01 0.15795E 06 0.25792E-03 0.40740E02 

0.48686E 04 0.39330E 00 : 0.95754E 00 -0.34560E-04 

— 0.90734E— 03 0.26868E 00 0.28180E 00 0.34869E 00 

0.50955E 00 0.42419E 00 0.81556E-02 0.46790E 00 


219 


0.59973E 00 0.0 0.33643E 00 C.33207E 00 

0.26048E 04 0.9291 IE 07 0.48936E-G4 0.79125E-Q6 

0.135I2E 01 0.L6110E 06 0.26451E-03 0.42613E 02 

G.64548E 04 0.42427E 00 0.11014E 01 G.91712E-05 

—0.944946—04 0.26872E 00 G.30629E 00 0.35087E 00 

0.62468E 00 0.455036 00 0.10041E-01 0.575926 00 

0.718936 00 0.37715E-01 0.34203E 00 0.32348E 00 

0.25879E 04 G. 933516 07 0.49517E-04 0.79354E-06 

0.13392E 01 0.161486 06 0.27736E-03 0.447S9E 02 

G.8Q5S66 04 0.45514E 00 0.126416 01 0.9G884E-05 

0.41139E— 04 0.26779E 00 -0.64252E-Q1 0.34154E 00 

0.739826 00 0.485816 CO 8.12218E-01 G.70O&9E 0© 

0-83814E 00 0*483166-01 0.33388E 00 C.30696E 00 

0.26I25E 04 0.92710E 0? 0.48672E-04 0.79019E-06 

0.13366E 01 0.16092E 06 O.29Q89E-03 0.46809E 02 

0.966006 04 0.485966 00 0.14407E 01 -0.16227 E-04 

-0. 147136-03 0.26772E 00 -0.989146-01 0.324446 00 

0.970076 00 0.54738E 00 0.17552E-01 9.100616 01 

0. 1-0 765 £. 01 0.0 0.332916 00 0.28306E 00 

0 .261556 04 D.92633E 07 0.485716-04 0.78979E-06 

0.135876 01 0.160856 06 0.30929E-03 0.497486 02 

0.« £28546 05 0.547666 00 0.18366E 01 -0.130366-04 

• 0.19794E-04 0.26804E 00 0.35037E-01 0.299216 00 

0®li632£ 01 C.59898E 00 0.23151E-01 0.132666 01 

0.1276SE 01 0.0 0.335426 00 0.27210E 00 

0.260786 04 0.92832E 07 0.48832E-04 0.79083E-06 

0.13533E 01 0.161036 06 0.32251E-03 0.519336 02 

0.155376 05 0.599426 00 0,221346 01 -0.224666-05 

0* 129346—03 0.267996 00 0.15065E-01 0.28754E 00 

0.I5495E 01 C.70190E 00 0.377876-01 0.21650E 01 

0. 167656 01 0.61 1266-02 0.33513E 00 0.252606 00 

0.260876 04 0.928096 07 0.488016-04 0.790716-06 

0.13540E 01 0. 16 10 IE 06 0.34897E-03 0.56186E 02 

0.20905E 05 0.70292E 00 0.30816E 01 -0.29927E-04 


220 


-0.74868E-04 0.26793E 00 -0.64Z66E-02 

0.19356E 01 0.8043IE 00 0.57674E-01 0, 

0.20765E 01 0.70041E-02 0.33286E 00 « 

0.26 15 6E 04 0.92629E 07 0.48565E-04 

0.13589E 01 0.16084E 06 0.37402E-Q3 

0.26269E 05 Q.80638E 00 0.41Q14E 01 

-0.223641-03 D.26793E00 -0.73719E-02 

0.23216E 01 0.90594E 00 0.83676E-01 0. 

0.24765E 01 0.320012-02 0.33153E 00 { 

0.26 198 E 04 0.92523E 07 0i48427E-04 

0.13617E 01 0.1 607 4E 06 0.397201-03 

- 0.31628E 05 0.90980E 00 " ' 0. 52.771 E 01 

— 0.94580E-04 0.26795E 00 -0^608632-03 

0.27074E 01 0.10064E 01 0.11664E 00 0. 

0.28765E 01 0.29892E-02 0.33086E 00 ( 

0.262172 04 0.92469E 07 0.48357E-04 < 

0 . 13632E 01 0.16069E 06 0.41910E-03 

0. 36986E 05 0.10132E 01 0.66108E 01 

— 0.59832E— 06 0.26796E 00 0. 1 76 7 5 E- 02 

0.30928E 01 0.11053101 0.15736E 00 0« 

0.32765E 01 0.617G3E-O2 0.32992E 00 C 

0.26246E 04 0.92394E 07 0.48259E-04 

0.13652E 01 0.16O62E 06 0.440212-03 

0.42341E 05 0.11164E 01 0.81025E 01 

-0. 1 1200E-03 C.26795E 00 0.86486E-03 

0.34778E 01 0.12020E01 0.20657E 00 0. 

0.36765E 01 0.10331E-01 0.32830E 00 ( 

0.26295E 04 G.92264E 07 0.48089E-04 

0.13688E 01 0.16050E 06 0.46079E-03 

0.47693E 05 0.12196E 01 0.97515E 01 

— Q.25064E— 03 0.26794E 00 • -0.17445E-03 

0.38622E 01 0.12958E 01 ; 0.26498E 00 0. 

0.40765E 01 0.146306-01 0.32605E 00 { 


0.26694E 00 

» 33043E 01 
5.23528E 00 
0. 78977 E— 06 
0.60158E 02 
— &«> 6 14QIE— 04 
0.2487IE 00 

.47922E 01 
5.22079E 00 
0.78922E-06 
0.638482 02 
-0.3Q533E-O4 
0.23344E 00 

. 66757E 01 
). 2091 0E 00 
0. 78893 E-06 
0.67346E 02 
0.25934E— 04 
v V 0.22 109E 00 

. 90002E 01 
). 19926E 00 
Q.78854E— 06 
0.70708E 02 
0- 42556E— 04 
0.210722 00 

. 1 L809E 02 
). 19044E 00 
0. 78 786 E— 06 
0.73956E 02 
0. 16471 E— 04 
0.201 44E 00 

. 15139E 02 
). 18238E 00 


221 


0.26364E 04 
0.13737E 01 
0.53040E 05 
— 0. 13630E— 03 


Q.92083E 07 
0.16032E 06 
0.13226E 01 
G.26794E 00 


0.47853E— 04 
0.48098E-03 
0.11557E 02 
— 0.406 79E— 03 


0-42459E 01 G.13859E 01 0.33318E 00 0- 

0.44765E 01 C.19464E-01 0.32324E 00 C 

0.26450E 04 0.91856E 07 0.47558E-04 

0.13799E 01 0. 160 10E 06 0.50098E-03 

0.58381E 05 0.14254E 01 0.13518E 02 

0.17279E-03 0.26795E 00 -0.15066E-03 


0.46287E 01 O.14715E01 0.41167E 00 0. 

0.48765E 01 0.25674E-01 0.31981E 00 ( 

0.26556c 04 0.91576E 07 0.47197E-04 

0.13875E 01 0.15982E 06 0.52091E-03 

0.63713E ®5 0.15280E 01 0.15630E 02 

0.36901E— 03 0-26795E 00 0.21338E-04 

0.50I03E 01 0.15514E 01 0.50083E 00 0. 

0.52765E 01 0.33494E-01 0-31565E 00 < 

0.26684E 04 0.91234E 07 G.46758E-04 

0.13968E 01 0.15946E 06 0.54088E-03 

0.69034E 05 0.16302E 01 C.17891E 02 

0.16770E— 03 0.26795E 00 0.89577E-04 

0.53906E 01 0.I6246E 01 0.60G86E 00 0. 

0.56765E 01 0.43072E-01 0.31066E 00 I 

0.26839E 04 0.90820E 07 0.46229E-04 

0.14081E 01 0.15902E 06 0.56103E-03 

0.74343E 05 0.17321E 01 0.20297E 02 

— G.4Q600E-03 0.26795E 00 0.T9145E-04 


0.57693E 01 
0.60765E 01 
0.2702 BE 04 
0.14217E 01 
0.79634E 05 
— 0.63293E— 03 


0.16898E 01 
0. 54459E- 01 
0.90324E 07 
0.15847E 06 
0-18336E 01 
0.26795E 00 


0.71182E 00 0. 

0.30477E 00 1 

0.45601E— 04 
0.581 59E- 03 
0.22843E 02 
0.27923E— 04 


0. 78691 E-06 
0.77113E 02 
0.31609E-05 
0.19297E 00 

19018E 02 
U17493E 00 
0.78571 E-06 
0.80207E 02 
0.2521OE-O4 
0.18516E 00 

23471E 02 
la 16803E 00 
0.78424E-06 
0.83249E 02 
0. 59871 E-04 
0.17793E 0© 

28514E 02 
)• 161 54E 00 
0.78244E— 06 
0.86249E 02 
0.6983OE— 04 
0.17116E 00 

. 34160E 02 
3.15537E 00 
0.78025E-06 
0.8921 BE 02 
0. 48495 E— 04 
0.16474E 00 

.40407E 02 
3. 1494 IE 00 
0.77762E— 06 
0.92162E 02 
0. 24378E— 04 
0.15855E 00 


222 


o • o • o ♦ 


61462E 0i 0» 1T458E OI 0.83353E 00 0.47239E 02 

• 64765E 01 0.67727E-O1 0.29792E 00 0.14357E 00 

Q. 272386 04 C.89740E 07 0.44867E-04 0. 7745IE-=-06 

0. 14377E 01 0. 15779E 06 0.6G279E-03 0.95113E 02 

G.84905E 05 0.193466 01 0.25522E 02 0. 11597E-04 

-0.39161E-G3 0.26794E 00 -Q.44366E-05 0.15250E 00 

652 10E 01 0.17913E 01 0.96561E 00 G.54634E 02 

. 66 765 E 01 G.83G80E— 01 0.29O10E 00 0.13781E 00 

0.27486E 04 G.89061E 07 0.44023E-04 0.77089E-06 

Q.14563E Oi 0.15697E 06 0.62479E-03 0.98071E 02 

0.90 152 E 05 0.203506 Oi 0-28327E 02 -0.64899F-C5 

— O.IO252E—03 0.26795E 00 -0.193756-04 O. 14655E 00 

58934E OI 0. 18251F Oi 0. 1IC75E 01 0.&2566E 02 

.72765c 01 0.1G067E 00 0.28129E 00 0.132G9E 00 

0. 27769c 04 0U88279E 07 O.43O63E-04 Q.76669E-06 

0.I4778E 0! 0. 15597E 06 G.64779E-03 O.L0104E03 

0.95368E 05 Q.2i348E01 C.31249E02 -0.28265E-04 

-0.168466-03 0.26795E 00 -0.19851E-04 0.14065E 00 

0.72633E 01 0.I8459E 01 0.12582E 01 0.70994E 02 

0.76765E OI 0.12Q59E DO 0.27150E 00 0.12638E 00 

Q-ESOSEE 04 0.87390E 07 0.41987E-04 0.76I9OE-O6 

G.I5023E Oi Oe 15479E 06 0.67202E-03 0.10402E 03 

0.10055E 06 0.22339E 01 0.34280E 02 -0.35497E-04 

— 0.29964E— 03 0.26795E 00 -0. 1 182GE-04 Q.13477E 00 

0.76304E 01 G.18525E 01 0.L4169E 01 0.79671E 02 

0.80765E 01 0.14301E 00 0.26077E 00 0.12Q64E 00 

0.28442E 04 0.86390E 07 0.40796E-04 0.75647E-06 

G.15300E 01 0.15338E 06 0.69772E-03 0.10702E 03 

0.105 68 E 06 0.23323E 01 0.3741 IE 02 -0.34790E-04 

-0. 15509 E— 03 0.26795E 00 -0.51113E-05 0.12S88E 00 

0.79945E OI 0.18441E 01 0.15822E 01 0.89152E 02 

0.84765E 01 0.16808E 00 0.24917E 00 0.11489E 00 

0.28831E 04 0.85274E 07 0.39491E-04 0.75038E-06 


223 


to S O 0 O 


0-1561 IE 01 0.15173E 06 0.72516E-03 0.11003E 03 

0. 11077 E 06 0.24298E 01 0.40630E 02 -0.22238E-04 

0.17993E-03 0.26795E 00 0.14718E-06 0.12295E 00 

83554E 01 0.18196E 01 0.17528E 01 0.98782E 02 

. 88765E 01 0.19593E 00 0.23679E 00 0-10909E 00 

0.29256E 04 G.84040E 07 0.38079E-04 0. 74360E-06 

0-15957E 01 0.14982E 06 0-75463E-03 0.11306E 03 

0.11580E 06 0.25265E 01 0.43928E 02 -0.769&2E-0S 

0.35769E— 03 0-26795E 00 0-44234E-05 0.11700S 00 

87131E 01 0.17786E 01 0-19271E 01 0.10871E 03 

-92765E 01 G.22682E 00 0-22373E 00 0.10327E GO 

0o297i4E 04 0.82689E 07 0.36567E-04 0.7361GE-06 

0-16338E 01 0-14761E 06 0.78646E-03 0.11609E 03 

0-12075E 06 0.26224E 01 0.47294E 02 0.143056-04 

G-32854E— 03 0-26795E 00 0.69551E-05 0.11101E 00 

90674E 01 0« 17204E 01 0.210336 01 0.11889E 03 

-96765E 01 0.26104E 00 0.21012E 00 0-97427E-01 

0.30205E 04 0.8122GE 07 0-34964E-04 0.72787E-06 

0.16758E 01 0.14509E 06 0.82101E-03 0.11912E 03 

0-I2563E 06 0.27173E 01 0.50719E 02 G.35702E-04 

0.172086-03 0.26795E 00 0.70083E-05 0-10499E GO 

0.94182E 01 G.16449E 01 0.22799E 01 0.12927E 03 

-0.10076E 02 0-29891E 00 0.19612E 00 0.91580E-01 

0.3Q725E 04 0.79635E 07 0.33283E-04 0.71890E-06 

0.17215E 01 0.14225E 06 0.85868E-03 0.12215E 03 

0.13042E 06 0.28113E 01 0.54190E 02 0.36361E-04 

—0- 10831E— 04 0.26795E 00 0.53534E-05 0.98955E-01 

0.97656E 01 0.15520E 01 0.24550E 01 0.13982E 03 

0.10476E 02 0.34079E 00 0.L8186E 00 0.85748E-01 

0.31273E 04 0.77937E 07 0.31536E-04 0.70918E-06 

0.17712E 01 0.139O7E 06 0.89996E-03 0.12515E 03 

0.13511E 06 0.29044E 01 0.57700E 02 0.21308E-04 

— 0.19843E-G3 0.26795E 00 0.33451E-05 0.92922E-01 


224 


o • o • O ® o * 


101 10E 02 0.14418E 01 0.26269E 01 0.15049E 03 

•10876E 02 0.38729E 00 0.16752E 00 0.79955E-01 

0.31846E 04 0.76129E 07 0.29739E-04 0.69870E-06 

0.18249E 01 0.13555E 06 0.94538E^03 0-12814E 03 

0.13969E 06 0.29966E 01 0.61239E 02 0-98072E-Q6 

-0.29816E-03 0.26795E 00 0.17802E-05 0-8691 3 E-01 

10450E 02 0.13149E 01 0.27939E 01 0.16126E 03 

• 11276E 02 0.43905E 00 0.15325E 00 0.74226E^0I 

0.32441E 04 G-74216E 07 0.27906E-04 0.68747E-06 

0-18828E 01 0.13169E 06 0.99557E-03 0.L3110E 03 

0.14415E 06 0.30878E 01 0.64797E 02 -0.15811E-04 

— 0.24717E— 03 0.26795E 00 0.83969E-06 0.80952E-01 

10787E 02 0.11716E 01 0.29543E 01 0.17210E 03 

-116766 02 0.49684E 00 0.13920E 00 0.68588E-01 

0.33055E 04 0-72205E 07 0.26054E-04 0.67549E-06 

0.19450E 01 0.12750E 06 0.10513E-02 0-13403E 03 

0.14847E 06 0.31782E 01 0-68367E 02 -0.23324E-04 

— 0- 19888E— 04 0.26795E00 -0.6564QE-07 0-75063E-01 

11121E 02 0.10129E 01 0.31067E 01 0-18299E 03 

-12076E 02 0.56207E 00 0-12552E 00 0.63070E-0I 

0-336856 04 0.70102E 07 0.24198E-04 0.66276E-06 

0-20116E 01 0. 12299E 06 0.11134E-02 0.13693E 03 

0-15264E 06 0-32677E 01 0-71942E 02 -0.159336-04 

0- 18133E— 03 0.26795E 00 -0.89380E-06 0.69276E-01 

0.11452E 02 0.83948E 00 0.32496E 01 0.19391E 03 

0-12476E 02 0.63588E 00 0-11233E 00 0.57698E-01 

0-343286 04 0-67914E 07 0-22354E-04 0.64929E-06 

0.20828E 01 0.11819E 06 0.11829E-02 0.13980E 03 

0-15666E 06 0.33563E 01 0-75515E 02 -0.21800E-05 

0.24267E-03 0-26795E 00 -0»16072E-05 0.636I7E-01 

0-11780E 02 0.65246E 00 0-33818E 01 0-20485E 03 

0.12876E 02 0-72004E 00 0.99760E-0I 0.52499E-01 

0-34982E 04 0.65649E 07 0-20537E-04 0.63510E-06 

0.21587E 01 0-11312E 06 0.12611E-02 0.14265E 03 


225 



0.16052E 06 0.34442E 01 0.79080E 02 0.19449E-04 

0.14947E-03 0.26795E 00 -0.20334E-05 0.58114E-01 

0.12105E 02 0-45292E 00 0.35021E 01 0.21581E 03 

0.13276E 02 0.81705E 00 0-87890E-0I 0.47497E-01 

0.35643E 04 0.63316E 07 0-18760E-04 0.62020E-06 

0.22397E 01 0.10781E06 0.13494E-02 0.14548E03 

0. 16420E 06 0.35313E 01 0.82631E 02 0.30970E-04 

— 0.28979E— 04 0-26795E 00 -0.2I943E-05 0.52792E-01 

0.12428E 02 0.24207E 00 0-36096E 01 0.22676E 03 

0.13676E 02 G-92886E 00 0.76801E-01 0.42714E-01 

Q.36308E 04 0.60922E07 0. 17037E-04 0.6Q460E-06 

0-23259E 01 0.10231E 06 0-14497E-02 Q-14832E 03 

0.16770E 06 0.36177E 01 0.86163E 02 0.20545E-04 

-0.28594E-03 0.26795E00 -0.21913E-05 0.47676E-01 

Q-12748E 02 G.21157E-01 0-37034E 01 0-23771E 03 

0.14076E 02 0.1O594E 01 0.66553E-01 0.38172E-0I 

0-36975E 04 Q.58480E 07 0.15381E-04 0.58833E-06 

0-24175E 01 0.96662E 05 0.15640E-02 0.15118E 03 

0- 171Q2E 06 0-37034E 01 0-89673E 02 0.64595E-06 

-0.44994E-03 0-26795E 00 -0.2Q440E-05 0.42791E-01 

0-13065E 02 -C.20852E 00 0-37828E 01 0.24866E 03 

0-14476E 02 0.L2125E 01 0.57180E-01 0.33890E-01 

0-37640E 04 0-55998E 07 0-13800E-04 0.57142E-06 

0.25150E 01 0» 90903E 05 0-16951E-02 0.15409E 03 

0.17415E 06 0.37885E 01 0.93155E 02 -0-23391 E-04 

— 0-40087E— 03 0-26795E 00 -0-18329E-05 0o38158E-01 

0- 13381E 02 — 0.44565E 00 0-38474E 01 0.25957E 03 

Oo 14876E 02 0.13913E 01 0-48707E-01 0-29882E-01 

0.38301E 04 0.53490E 07 0-12306E-04 0-55392EH56 

0-26184E 01 0.85093E 05 0-18460E-02 0.15708E 03 

0.17708E 06 0.38731E 01 0o96607E 02 -0.49206E-04 

— 0-13600E-03 0.26795E 00 -0-13312E-05 0.33796E-01 

0.13694E 02 -0.68889E 00 0.38967E 01 0-27045E 03 
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0»15276E 02 0.16026E 01 0«>4il38E^Ol 0.26166E-01 

0.38953E 04 0.50970E 07 0.10908E-04 0. 53588 E-06 

0o27281E 01 0.79289E 05 0®20202E-02 0-I6018E 03 

0.17982E 06 0.39571E 01 0.10003E 03 -0.59158E-04 

0.13195E-03 0. 26795E 00 -0«. 57991 E-06 0-29727E-01 

14006E 02 — 0.93694E 00 0.39306E 01 0.28127E 03 

• 15676E 02 0.18499E 01 0.34457E-01 0.22751E-01 

0.39594E 04 0.48453E 07 0.9611 0E-05 0.51739E— 06 

0.28440E 01 0.73549E 05 Q.22220E-02 0.16343E 03 

0.18237E 06 0.40407E 01 0-10341E 03 -0.42818E-04 

0.34986E-03 0.26795E 00 0.28027E-06 0.25965E-01 

14317E 02 -0.11885E 01 0.39489E 01 0.29201E 03 

-16076E 02 0.21343E 01 0.28645E-01 0.19648E-01 

0.40218E 04 0.45962E 07 0.84229E-05 0.49857E-06 

0.29661E 01 0.67944E 05 0.24558E-02 0.16686E 03 

0o 18472E 06 0.41239E 01 0.10675E 03 0.95799E-06 

0.45003E— 03 G.26795E 00 0.88889E-06 0.22525E-01 

14626E 02 -0.14422E 01 0.39518E 01 0.30262E 03 

* 16476E 02 0.24572E 01 0.23663E-01 0.16866E-01 

0o40821E 04 0.43520E 07 0.73484E-05 0.47961E-06 

0.30939E 01 0.62544E 05 0.27259E-02 0.17049E 03 

0« 18690E 06 0o42067E 01 0.11005E 03 0.22452E-04 

0.32487E— 03 0.26795E00 0.11548E-05 0.19421E-01 

0» 14934E 02 -0.16969E 01 0.39394E 01 0.31304E 03 

0.16876E 02 0.28032E 01 0.19464E-01 0.14407E-0I 

0.41395E 04 0.41158E 07 0.63914E-05 0.46072E-06 

0.32262E 01 0.57426E 05 0.30357E-02 0.17433E 03 

O.iaaOOE 06 0.42893E 01 0.1L330E 03 0.22726E-04 

0.46111E— 05 0.26795E 00 0.11673E-05 0.1666IE-0I 

0.15242E 02 -0.19513E 01 0.39120E 01 0.32318E 03 

0.17276E 02 0.31402E 01 0.15999E-0I 0.12275E-0I 

0.41933E 04 0.38916E 07 0.55562E-05 0.44227E-06 

0.33610E 01 0.52680E 05 0.33853E-02 0.17834E 03 

0.19073E 06 0.43717E 01 0.11651E03 0.83962E-05 
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-0.29Q77E-03 


0.26795E 00 


0.11453E-G5 


0.14254E-01 



0.15549E 02 -0.22043E 01 0.38702E 01 0.33294E 03 

0* 176T6E 02 0o34162E 01 0.13205E-01 G.10469E-0S 

C.42426E >34 0.36839E 07 0.48443E-05 O.42469E-06 

G.34950E 01 0.48394E 05 0.37698E-02 0.18244E 03 

G.1924IE 06 0.44539E 01 0.11966E 03 -O.59017E-G5 

-0.30794E-03 0-26795E 00 0.T1174E-05 0.12203E-SS. 

G.15855E 02 -0.245.46E 01 0.38I47E 01 0.'3*219E 03 

0.18076E 02 0.35368E 01 0.11016E-01 0.89796E-02 

G.42862E 04 C..34960E 07 0.42562E-05 0.40853E-06 

0.36235E 01 0.44655E 05 Q.41754E-02 0.18645E 03 

0. 19396E 06 0.45362E 01 0.12276E 03 -G.23416E-G4 

— G.JL0769E— 03 0.26795E 00 0.12089E-05 0.10503E-01 

0.L6163E-02 -0.27013E 01 0.37462E 01 0.35074E 03 

0.I8476E 02 0.33855E 01 0.93693E-02 0.77918E-02 

Q.43229E 04 0-33399E 07 0.37913E-05 0.39443E-06 

G.374G1E 01 0.41552E 05 0.45755E-02 0.19012E 03 

0.I9540E m 0.46185E 01 0.12579E 03 -0.92975E-05 

G.13926E-03 0.26795E 00 0.12131E-05 0.91403E-02 

0* £647 IE 02 -0.29433E 01 0.36657E 01 0-35840E 03 

G.18876E 02 C.28747E 01 0.81995E-02 Q.68808E-02 

0.43517E 04 0.32150E 07 0.34468E-05 G.38308E-06 

G.38374E 01 0.39155E 05 0.49325E-02 0-19313E 03 

0-19674E '06 0.47011E 01 0.12876E 03 0.19789E-05 

0.26183c- 03 0.26795E 00 0.10333E-05 0.809D3E-02 

0.16780E 02 — 0.3 1796E 01 0.35744E 01 0.36501E 03 

0.192T6E 02 0.20 l 27E 01 0.74435E-02 0.62034E-02 

0.43718E 04 C.31274E 07 0.32167E-05 0. 37497 E-06 

0.39088E 01 0.37503E 05 0.52044E-02 8.19518E 03 

0.19802E 06 0.47840E 01 0.I3165E 03 0. 11248E-04 

0.28066E— 03 0.26795E 00 0.71640E-06 0.73057E-02 

0.I6936E 02 -0.32955E 01 0.35251E 01 0.36768E 03 

0.19476E 02 0.14950E 01 0.72016E-02 0.59313E-02 
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0.43785E 04 0.30980E 07 0.31417E-05 0.37223E-06 

0.39333E 01 0.36955E 05 0-52978E-02 0ol9578E 03 

0. 19864E 06 0..48256E 01 0.13307E 03 0.D0143E-04 

0.25678E— 03 0.26795E 00 0.58670E-06 0.69891E-02 

0o I7248E 02 -0.35221E 01 0..34201E 01 0.37278E 03 

0.19876E 02 0.59421E 00 0.69358E-02 0o54916E-02 

0..43861E 04 0o30649E 07 0o30584E-05 0.36913E-06 

0.39613E 01 G.36340E 05 0„53987E-02 0<>19619E 03 

0«> 19986E 06 0.49094E 01 0.13584E 03 0.49655E-05 

0. 15818E— 03 0c 26795E 00 0.53421E-06 0.64750E-02 

0.17563E 02 -0.37415E 01 0.33074E 01 0.37684E 03 

9.20276E 02 0.48018E-01 0«>68631E-02 0.51098E-02 

0.43882E 04 G.30557E07 0.30355E-05 0.36827E-06 

0.39692E 01 0.36170E 05 0.54153E-02 0ol9587E 03 

0. 20106E 06 0.49937E 01 0.13852E 03 ~ 0 « 15922E-05 

0.57871E-04 0.26795E 00 0.64454E-06 G.60259E-02 

0o 17881E 02 — 0- 39535E 01 0-31881E 01 0.38046E 03 

0.20676E 02 0.0 0.68804E-02 G.50406E-02 

0o4387?E 04 0.30579E 07 0.30409E-05 0. 36847E-06 

0.39673E 01 0.36210E 05 0.53918E-02 0ol9524E 03 

0.2Q227E 06 0.50788E 01 0.14U2E 03 -0.22789E-05 

0..37167E— 04 0.26795E 00 0.76284E-06 0o59440E-02 


17960E 02 — 0.40053E 01 0.31574E 01 0.38135E 03 

« 20776E 02 0.0 0.68883E-02 0«50492E-02 

0.43874E 04 0.30589E 07 0o30434E-05 0.36857E-06 

0.39664E 01 0.36229E 05 0.53845E-02 0.19507E 03 

0.20257E 06 0.51002E 01 0.14175E 03 0.10371E-04 

0.52652E— 04 0.26795E 00 0.78886E-06 0.59540E-02 


O^lSOOOE'vK? — 0.40311E 01 0.31418E 01 0<,38179E 03 

0.20826E 02-^ 0.0 0.68919E-02 0.50532E-02 


0.43873E 04 *w?0593E 07 0.30446E-05 0.3686IE-06 

0.39660E 01 0.3(3^13,2^ 05 0.53810E^02 0.19499E 03 

0.20272E 06 0 . 5 1 1 0.14207E 03 0.12409E-04 

0.64107E— 04 0.26795E Off-*.,. 0.7697 IE-06 0.59587E-02 

X 



o • o « 


0-18080E 02 — G.40823E 01 0.31106E 01 G.38268E 03 

G.2Q926E 02 0.0 0.68981 E-02 G.50606E-02 

0 .4367 IE 04 0.30601E 07 0.30465E-05 Q.36868E-06 

0-33654E 01 0.36252E 05 0.53745E-G2 0. 19484E 03 

0.20302E 06 0.51323E 01 0.14269E 03 0.16032E-04 

0.92141E-04 0.26795E 00 0.78392E-06 0.59674E-02 

* 18241E 02 -G.41832E 01 0.30470E 01 0.38449E 03 

. 21 126E 02 0.0 0. 69062 E-02 G.50726E-02 

0.43869E 04 0.30611E 07 0.30491E-05 C.36878E-0S 

G.39645E 01 0.36271E 05 0.53637E-02 0.19455E 03 

0.20363E 06 0.51753E 01 0-14393E 03 0.20887E-04 

0.15526E— 03 0.26795E 00 0.72558E-06 Q.59813E-02 

18563E 02 — 0.43799E 01 0.29160E 01 0.38818E 03 

.21526E 02 0 . 1 1134E-01 0.69096E-02 G.51242E-02 

0.43868E 04 0.3O616E 07 0.30501E-05 0.36882E-O6 

0.3964IE 01 0.36279E 05 0.53483E-02 0.194Q3E 03 

0.20484E 06 0.52618E 01 0.14635E 03 0.28419E-04 

0.25596E— 03 0.26795E 00 0.50975E-06 0.60422E-02 

0» 18889E 02 -C.45695E 01 0.27804E 01 0.39I93E 03 

0.21926E 02 0.33701E-01 0.69032E-02 0.51728E-02 

0. 43870 E 04 0.30608E 07 0.30481E-05 0.36874E-06 

0.'39648E-0I Q.36264E 05 0.53379E-02 0-19357E 03 

0.20605E 06 0.53489E 01 0.14868E 03 0.35014E-04 

0.24351E— 03 0.26795E 00 0.24515E-06 0.60995E-02 

0.19216E 02 — 0.47522E 01 0.26405E 01 0.39570E 03 

0.22326E 02 Q.30IUE-01 0.68944E-02 0.51739E-02 

0.43872E 04 Q.30596E 07 0.30454E-05 0.36864E-06 

0.39658E 01 0.36243E 05 0.53291E-02 0.19314E 03 

0.20723E 08 0.54365E 01 0.15094E 03 0.34164E-04 

0.6626 3E— 04 0.26795E 00 -0.38373E-0? 0.61010E-02 

0.19380E 02 -0.48411E 01 0.25692E 01 0.39757E 03 

0.22526E 02 0.20653E-01 0.68911E-02 0.51619E-02 

0.43873E 04 0-30592E 07 0.30443E-05 0.36860E-06 
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0.39661E 01 0.362366 05 0.53244E— 02 0.19293E 03 

0»20786E 06 0.54806E 01 0.15204E 03 0.30874E-04 

— 0.79056E— 04 0.26795E 00 -0.16857E-06 0.60869E-02 

0. 19710E 02 -0.50139E 01 0.24239E 01 0.40131E 03 

0. 22926E 02 0.50047E-02 0.68879E-02 0.51371E-02 

0.43874E 04 0.30588E 07 0.30433E-05 0.36856E-06 

0.39665E 01 , 0.36228E 05 0.53142E-02 0.19252E 03 

0o 20907E 06 0.55690E 01 0.15420E 03 0.39408E-04 

— 0.29880E— 03 0.26795E 00 -0.53514E-06 0.60577E-02 

0-20042E 02 -0.51803E 01 0.22755E 01 0.40502E03 

0.23326E 02 0.0 0.68879E-02 0.51265E-02 

0.43874E 04 0.30588E 07 0.30433E-05 0.36856E-06 

0.39665E 01 0.36228E 05 0.53032E-02 0.19212E 03 

0.21027E 06 0.56580E 01 0.15629E 03 0.328786-04 

— 0.37869E— 03 0.26795E- 00 -0. 10113E-05 0.6045 IE -02 

0.20376E 02 -0.53405E 01 0.21244E 01 0.4087 IE 03 

0.23726E 02 0.0 0.68896E-02 0.51375E-02 

0.43874E 04 0.30590E 07 0.30438E-05 Q.36858E-06 

0.39663E 01 0.36232E 05 0.52922E— 02 0.19175E 03 

0. 21148E 06 0.57475E 01 0. 15831E 03 0.19463E-04 

— 0.33067E— 03 0.26795E 00 -0.14882E-Q5 0.60582E-02 

" 20712E 02 — 0.54947E 01 0. 19708E 01 0.41241E 03 

. 24126E 02 0.0 0.6891 7E-02 0.51483E-02 

0.43873E 04 0.30593E 07 0.30445E-05 0.36861E-06 

0.39661E 01 0.36237E 05 0.52816E-02 0.19139E 03 

0.21269E 06 0.58374E 01 , 0.16027E 03 0.17064E-05 

— 0.19673E-03 0.26795E 00 -0.17563E-05 0.60708E-02 

0.21049E 02 — 0.56430E 01 0.18153E 01 0.41611E 03 

0.24526E 02 0.0 0.68931E-02 0.51581E-02 

0.43873E 04 0.30595E 07 0. 30449E-05 0. 36862 E-06 

0.39659E 01 0.36240E 05 0.52719E-02 0.19105E 03 

0.21390E 06 0.59278E 01 0.16217E 03 -0.986466-05 

-0.54979E-04 0.26795E 00 -0. 16893E-05 0.60823E-02 
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0.21388E 02 -0 

0.24926E 02 
0.43873E 04 
0.39659E 01 
0-21510E 06 
0.36107E— 05 


•57858E 01 
0.0 

0-30595E 07 
0.3624LE 05 
0.60186E 01 
0.26795E 00 


0.16578E 01 
0.68936E-02 
0. 30451E— 05 
0. 52631E-G2 
0.16401E 03 
-0.11835E-05 


0o 4 198 IE 03 
0.51668E-02 
0. 36863E— 06 
0-19074E 03 
— 0.42362E— 05 
0.609 26E-Q2 


0.2I728E 02 -0.59231E 01 0.14989E 01 0.42352E 03 

0. 25325E 02 G.97009E-03 0.68935E-02 0.51817E-02 

0.43873E 04 0.30595E 07 0.30451E-05 0.36863E-06 

0.39659E 01 0.36241E 05 0.52550E-02 0.19045E 03 

G.21631E 06 0.61098E 01 0.16580E 03 -0.11520E-04 

—0.9076 IE-05 0.26795E 00 -0.38535E-06 0.61102E-02 

0.22070E 02 -0.60552E 01 0.13386E 01 0.42723E 03 

0.25725E 02 0.16920E-02 0.68931E-02 0.51921E-Q2 

0.438736 04 0.30595E 07 0.30449E-05 0.36862E-06 

0.39659E 01 0.36240E 05 0.52477E-02 0.190I8E 03 

0.21752E 06 0.62014E 01 0.16753E 03 0.64796E-05 

— 0.65084E— 04 0.26795E 00 0.50460E-06 0.61225E-02 


*** CALCULAT I GN5 FOR BETA =0.10000E 01 HAS TERMINATED 
AFTER 93 INTEGRATIONS ** 


FOR STREAMLINE NO. 2 WHERE BETA = 0.45000E 02 DEG. 
THE CALCULATIONS ARE AS FOLLOWS: 


X 

S 

UE 

MACH NO 
RN1 

D2R/DPHI2 

0. 33368 E-01 
0.41158E— 01 


Y 

BBAR 

HE 

RN/FT 

R 

DR/DX 

0.15237E 00 
0.72033E 00 


Z 

P/PS 
RHOE 
MOM THK 
PHI 

D2R/DX2 

0.26340E-0! 
0.99219E 00 


H 

QW/QWS 

HUE 

RNM 

DR/DPHI 

QB/QSB 

0.37250E— 01 
0.89825E 00 


o o o 0 O 0 


0®23824E 03 0.12655E 08 0®10596E- 03 0®95225E-06 

0ol0589£ 00 0®26509E 05 0®15225E-03 Q*4G361E 01 

0.93798E-02 0.15463E 00 0®98078E 01 -0.46050E-05 

0.45506E-04 0®22093E 01 -0®38022E 02 0®89865E 00 

3 4631 E— 01 0® 15470E 00 0®28954E-01 0®40954E-01 

° 44883E-01 0®65260E 00 0.990I1E 00 0.92819E 00 

0o26816E 03 0.12648E 08 0ol0580E-03 0®95192E-06 

0® 11922E 00 Qo29804E 05 0®14656E-03 0®43682E 01 

0®87518E 01 0®15739E 00 O 0 IO 6 OIE 02 0.53337E— 05 

0®67907E-G4 0®21625E 01 -0o36060E 02 0o92872E 00 

37265E— 01 0® 15935E 00 0«34150E-01 0®48269E-01 

o 52333E-01 0® 58865E 00 0®98537E 00 0®95743E 00 

0®32647E 03 0®12630E 08 0. 1Q544E-03 0<- 95 1 1 5E— 06 

0ol4525E 00 0o36190E 05 0®1412QE-Q3 0®51101E 01 

0o29249E 02 0®16296E 00 0®12096E 02 0®19980E-04 

0.1Q478E— 03 0»20724E 01 -0®32492E 02 0-95824E 00 

42947E-01 0® 16849E 00 0®44452E-0I 0®62665E-0! 

*> 67233E— 01 0-54975E 00 0®97356E 00 0®96995E 00 

0®43988E 03 0«12587E 08 0®10453E-Q3 0®94923E — 06 

0ol9604E 00 0®48441E 05 0®13841E-03 0»67048E 01 

0®81850E 02 0.17425E 00 0®14779E 02 0.32529E-04 

0o 14706E-03 Go 19056E 01 -0® 26595E 02 0.97144E 00 

0o49162E— 01 0o 17738E 00 0.54661E-01 0®76856£=01 

0o 82133E— 01 Oo 54158E 00 0.95870E 00 0.96478E 00 

0o55125E 03 0.12532E 08 0.10339E-03 0®94678E-06 

0o24622E 00 0«60198E 05 0®13836E-03 0®83287E 01 

0® 14934E 03 0®18561EOO 0®17127E 02 0.28250E-04 

0® 14610E-03 0o 17555E 01 -0®22023E 02 0®96710E 00 

0®55892E-01 0® 18601E 00 0®64778E-01 0®90875E-0l 

0-97033E-01 0®54277E 00 0„94089E 00 0.95382E 00 

O 066 I 68 E 03 0®12465E 08 0®10202E^03 0„94381E-06 

0® 29633E 00 0®71520E 05 0®13909E-03 0®99478E 01 

0®23116E 03 0®19696E 00 0.19201E 02 0.13069E-04 

0®98928E— 04 0®16202E 01 -0®18452E 02 0®95713E 00 
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0.70859E-0I 0 

0.12683E 00 
0.88168E 03 
Q.39758E 00 
0.43562E 03 
-0.84918E-04 


.20237E 00 
0.55873E 00 
0.12295E 08 
0.92833E 05 


0.84673E-01 
0.89679E 00 
0.98577E-04 
0.14167E-03 


0-11835E 00 
0.92289E 00 
Q.93623E— 06 
0.13 152E 02 


0.219376 00 
0.13866E 01 


G.22705E 02 
-0.13395E 02 


—0. 54221E-05 
0.92859E 00 


0.87755E— 01 C.21748E 00 0. 10400E 00 0* 14487E 00 

0.15663E 00 G.57968E 00 0.84245E 00 0.88160E 00 

0.11012E 04 0.12077E 08 0.94272E-04 0.92643E-06 

0.5G1G2E 00 0.11206E 06 0.14551E-03 0.163G6E 02 

0.69052E 03 0.24107E 00 0.25558E 02 -0.12364E-04 

— 0.20395E— 03 0.11917E01 -0.10108E02 O.89O12EG0 


0.10647E 00 0.23124E 00 0.12266E 00 0.17027E 00 

0.18643E 00 C.60494E 00 0.78021E 00 0.83092E 00 

0.L3178E 04 0.11815E 08 0.89243E-04 O.91451E-06 

0.60616E 00 0.12860E 06 0.15083E-03 0.19396E 02 

0.98994E 03 0.26176E 00 0.27943E 02 -0.26751E-04 

— 0-75924E— 04 0.10259E 01 -0.78617E 01 0.84247E 00 


0.12688E 00 0.24353E 00 0.14055E 00 0.19442E 00 

0.21623E 00 0.63700E 00 0.71265E 00 0.77344E 00 

0.15298E 04 0.11514E 08 0.83653E-04 0.90059E-06 

0.71264E 00 0.1421OE 06 0.15768E-03 0.22406E 02 

0.1 3267 E 04 0.28118E 00 0.29990E 02 -0.42208E-G4 

0. 2162 IE— 03 0.88226E 00 -0.62505E 01 0.78800E 00 


0.13768E GO 0.24909E 00 0.14918E 00 0.20599E 00 

0.23113E 00 0.65595E 00 0.67767E 00 0.74288E 00 

0. 16337 E 04 0.11349E 08 0.80698E-04 0.892926-06 

0.76677E 00 0.14765E 06 0.16172E-03 0.23877E 02 

0.15067E 04 0.29035E 00 0.30917E 02 -0.40671E-04 

0.29999E— 03 0.81711E 00 -0.56040E 01 0.75887E 00 


0.16037E 00 
0.26093E 00 
0.18399E 04 


0.25899E 00 
0.72884E 00 
0.109916 08 


0.16575E 00 
0.60573E 00 
0.74482E— 04 


0.2281 IE 00 
0.68070E 00 
0* 87600E— 06 
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0o87750£ 00 0ol5644E 06 0®m02E-03 0o26754E02 

0® 18849E 04 0.30749E 00 0.32619E 02 -0.13269E-04 

0e 21497E-03 0®69779E 00 -0.46729E 01 0.69939E 00 

0.18440E 00 0.26723E 00 0®18132E 00 0.24887E 00 

0.29073E 00 0.87197E 00 0.52873E 00 0®61672E 00 

0.20549E 04 0.10572E 08 0®67589E-04 0®85583E-06 

0o99928E 00 0.16229E 06 0.18226E-03 0.29579E 02 

0.22812E 04 0.32294E 00 0.34158E 02 0. 15782E-04 

— 0o20934E— 04 0.59040E 00 -0.42176E 01 0.63798E 00 

0® 19687E 00 0.27072E 00 . 0.18868E 00 0.25874E 00 

0® 30563E 00 0.97949E 00 0.48744E 00 0-58282E 00 

0.21699E 04 0®10329E 08 0®63775E-04 0.84394E-06 

0® 10675E 01 0.16397E 06 0.1890IE-03 0®3Q993E 02 

0-24839E 04 0«32999E 00 0.34875E 02 0-19524E-04 

— 0®5373QE— 04 0.54031E 00 -0®41455E 01 0«>60532E 00 

0.20963E 00 0.27380E00 0.19574E 00 0.26824E00 

0.32053E 00 0. I1184E 01 0.44389E 00 0®54666E 00 

0.22920E 04 Q.10057E 08 0®59652E-04 0.83040E-06 

0®11428E 01 0.16465E 06 0®19695E-03 0®32428E 02 

0o26880E 04 0.33657E 00 0o3556lE 02 0.17536E-04 

0.61657E— 04 0®49221E 00 -0.41554E 01 0.57032E 00 

0o23592E 00 Q.27871E00 0.20886E00 0o28598E00 

0o35033E 00 0«13100E 01 0.35555E 00 0.46463E00 

0®25474E 04 0®94391E 07 0®50908E-04 0.79895E-06 

0.L3110E 01 0. 16231E 06 0-21784E-03 0®35359E 02 

0o30950E 04 0.34828E 00 0.36847E 02 0® 12161E-04 

0o64460E— 03 0-39991E 00 -0.37209E 01 0®48972E 00 

0«26312E 00 0.28187E 00 0®2Z060E 00 0®30159E 00 

0.38G13E 00 0.11829E 01 0.28348E 00 0.38350E 00 

0.27699E 04 0.88475E 07 0®43302E-04 0.76775E-06 

0.14724E 01 0.15623E 06 0«24439E-03 0®38180E 02 

0®34910E 04 0.35794E 00 0.38048E 02 0.17433E-04 

0.73697E-03 0.31124E 00 -0.24684E 01 0o40823E 00 
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0.27003E 00 0.28237E 00 0.22333E 00 0.305116 00 

0.38758E 00 0.10566E Oi 0.27029E 00 0.365626 00 

0.28127E 0% 0.87279E 07 0.41854E-04 0.76130E-G6 

0.15054E Ci 0.15463E 06 0.251026-03 0.38816E 02 

0.35S75E 04 0.36001E 00 C.38340E 02 0.13447E-04 

G.S0S45E-03 0.28948E 00 -0.20232E 01 0.38998E 00 

0.2735GE 00 0*282586 00 0.22466E 00 0.3068IE 00 

0-39131E 00 0.97446E 00 C.26461E 00 0.357186 CO 

0.28314E 04 0.86751E 07 0.41223E-04 Q.75844E-06 

0*152006 Cl 0.1539QE 06 0.254176-03 0.39117E 02 

0.363536 04 0.361006 00 0.38486E 02 0.97881E-05 

0.34126E— 03 0.27865E 00 -0.17804E 01 C.38I31E 00 

0.27S24E DO (.282676 00 0.22532E 00 0-307646 00 

O. 393176 GO 0.930196 00 C.26201E 00 0.3531 IE 00 

0.2 840 IE 04 0.865066 07 0.409336-04 0.757116-06 

G. 152686 01 0.15355E 06 0.25570E-03 0.39262E 02 

0.365926 04 0. 361486 00 0.38558E 02 0. 768366-05 

G. 242456—03 0.273486 00 -0.16583E 01 0.37712E 00 

0.276736 OO 0.282836 00 0.22662E 00 0.309296 00 

0.396906 00 O. 859706 00 0.257236 00 0.34571E GO 

0*285.606 04 0.860536 07 0. 403996-04 0.754646-06 

0.153946 01 0.152896 06 0.25864E-03 C.39544E 02 

0.370686 04 0.36242E 00 0.38704E G2 0.39408E-06 

0.228056-04 C-26677E 00 -0.14708E 01 0.36951E 00 

0.28047S 00 0*282906 00 0.227276 00 0.31010E 00 

0.398766 00. C. 834936 00 0.255016 00 0.342396 00 

0.2S634E 04 0.85840E 07 0.40150E-04 0.75348E-06 

0. I5453E 01 0.15258E 06 0.26006E-03 0.39680E 02 

0.373056 04 0.362896 00 G.38776E 02 -0. 313886-05 

-0.920926-04 0.265346 00 -0.14070E 01 C„ 366096 00 

0.283966 00 0.28305E 00 0.22857E 00 0.3I173E 00 

0.402486 00 C.80U3E 00 0.25081E 00 0.336336 00 

0.287766 04 C.85433E 07 0.39676E-04 0.75126E-06 

D.15566E 01 0.151976 06 G. 262846-03 0.399456 02 
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0. 37777 E 04 0.36381600 0.38922E 02 -0.1U76E-04 

— 0.31796E— 03 0.26519E 00 -0.13226E 01 0.35986E 00 

0.29094E 00 0.28332E 00 0.23117E 00 0.31498E 00 

0.40993 E 00 0.73537E00 0.24320E 00 0.32515E00 

0.29035E 04 0.84684E 07 0.38812E-04 0.74715E-06 

0.15776E 01 0.15083E 06 0.26817E-03 0.40448E 02 

0.38717E 04 0.36566E00 0.39212E02 -0.29603E-04 

— 0.73510E-03 0.26551E 00 -0.U692E 01 0.34833E 00 

30488E 00 0.28380E 00. 0.23641E 00 0.32148E 00 

• 42483E 00 0.59900E 00 0.23083E 00 0.30592E 00 

0.29464E 04 0.83431E 07 0.37392E-04 0.74022E-06 

0. 161 29E 01 0.14883E 06 0.27780E-03 0.41345E 02 

0.40577E 04 0.36937E 00 0.39795E 02 -0. 66938 E-04 

-G.14197E-02 0.266O9E OO -0.88499E 00 0.32843E 00 

33270E 00 0.28448E GO 0.24706E 00 0. 33455E 00 

. 45463E 00 C.33C93E 00 0.21537E 00 0.27755E 00 

0.30014E 04 0.81795E 07 0.35586E-04 0.73110E-06 

0.16593E Gl €■« 14&09E 06 0.29269E-03 0.42760E 02 

0.44236E 04 0.37678E 00 0.40973E 02 -0.13529E-03 

— 0.21677E-02 0.26709E 00 -0.40455E 00 0.29880E 00 

8806E 00 0.28445E 00 0.26914E 00 0.36104E 00 

. 51423E 00 0.0 G.208G7E 00 0.24556E 00 

0.30280E 04 0.809936 07 0.34720E-04 0.72659E-06 

0. I6823E 01 0.14469E 06 0.30780E-03 0.44535E 02 

0.51439E 04 0.39160E 00 0.43416E 02 -0.18554E— 03 

— 0.16455E— 02 0.26849E 00 0.22982E 00 0.26472E 00 

0.49718E 00 0.27753E 00 0.31650E 00 0.41590E 00 

0.63 344 E 00 0.0 0.21810E 00 0.23923E 00 

0.29915E 04 0.82089E 07 0.35907E-04 0.73275E-06 

0.16509E 01 0.14659E 06 0.32059E-03 0.46996E 02 

0.65916E 04 0.42095E 00 0.48753E 02 — 0.191 77E— 04 

-0.25902E-03 0.26896EOO 0.36986E00 0.25742E 00 

0.524116 00 0.27424E 00 G.32883E 00 0.43014E 00 
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0.663246 00 Q.94247E-01 0.21691E 00 0.24175E 00 

0-29958E 04 G.61961E 07 0.35767E-04 0.73203E-06 

0.16546601 0.14638E 06 0.32598E-03 0.47715E02 

G.69555E 04 0.42818E00 0.50172E 02 0.110136-04 

-0.64561E-03 0.26863E 00 0.23983E 00 0.26019E 00 

Q.57757E 00 0.26583E 00 0.35379E 00 D.45929E 00 

0-72284E 00 0.25041E 00 Q.20846E 00 G.23412E 00 

G.30265E 04 O.S1036E 07 0.34767E-04 0.72684E-06 

0.168106 01 0. 14477E 06 0.34038E-G3 0.49275E 02 

0.76790E 04 0.44253E 00 C.53079E 02 0.12524E-04 

— 0* i0433£-G2 0.26807E 00 0.422466-01 G.2S237E 00 

0.63055E 00 0.25505E 00 0.37887E 00 Q.48909E 00 

0.78244E 00 0.36544E 00 G.19509E 00 0.22066E 00 

0.30764E 04 0.795I6E 07 0.33159E-04 G.71822E-06 

0.17250E 01 0.14203E 06 0.35830E-03 0.50889E 02 

0*839166 ©4 0.45672E 00 0.56052E 02 -G.18808E-04 

-0. 640746-03 0.26771E 00 -G.76486E-01 0.23848E 00 

0.683G7E 00 0.24 1986 00 0.40382E 00 9.51924E 00 

G.S4204E 00 0.44324E 00 0.17972E OG 0.20501E 00 

0.31357E 04 Q.77673E 07 0.31271E-04 0.70766E-06 

0. 17? 90S 01 0 . 13856E 06 0.37875E-Q3 0.52480E 02 

0.90886E 04 0.47077E 00 0.590696 02 -Q.16674E-04 

—0. I2343E— 03 0.26758E 00 -0.12266E 00 0-22227E 00 

0.73514E 00 0.22667E 00 0.42844E 00 0.54937E 00 

0.9C164E 00 0.48295E 00 G.16444E 00 0.18915E 00 

0.31972E C4 0.75726E07 0.29348E-04 0.69635E-06 

0. 1837 OE 01 0.13475E 06 0.40078E-03 0.54003E 02 

0.97674E 04 0.48471E00 0.62118E02 0.11722E-04 

— 0.582O4E— 04 0.26766E 00 -0.10442E 00 0.20576E 00 

0.78677E 00 C.20918E 00 0.45253E 00 0.57913E 00 

0.96124E 00 0.49766E 00 G.15053E 00 0.17443E 00 

0.32557E 04 Q.73838E 07 0.27552E-04 0.68523E-06 

0.18944E 01 0. 13091E 06 0.42346E-03 0.55434E 02 

0.10427E 05 0.49853E 00 0.65192E 02 0.67213E-G5 


233 


9 O 


-0.11233E-03 0®26782E 00 -0®48663E-01 0.19036E 00 

0o88867E 00 0.16774E 00 0.49836E 00 0.63705E 00 

0. 10804E 0 1 0.56771E 00 0»12628E00 0.14928E 00 

0o33649E 04 0-70224E 07 0®24304E-04 0.66350E-06 

0.20077E 01 0. 12325E 06 0.47148E-03 0-5811 IE 02 

0.1169QE 05 0.52583E 00 0.71398E02 -0.17080E-04 

0.18573E— 03 0-26799E 00 0.19311E-01 0.16393E 00 

0-98877E 00 Q.11816E 00 0-53988E 00 Q.69247E 00 

0-11996E 01 0.70654E 00 0.10464E 00 0.12725E 00 

0.34723E 04 0-66552E07 0.21250E-04 0.64079E-06 

0-21282E 01 0.11515E 06 0.52745E-03 0.60736E 02 

0- 12874E 05 0.55266E 00 0.77654E 02 0.31395E-04 

0.88272E-04 0.26804E 00 0. 34156E-01 0.14063E 00 

0® 10872E 01 G-61346E— 01 0.57578E 00 0.74531E 00 

0.13188E 01 0.92244E 00 0.84457E-01 Q.10658E 00 

0® 35843E 04 0.62599E 07 0«18234E-04 0.61556E-06 

0-22651E 01 O.10617E 06 0.59776E-03 0.63467E 02 

Do L3974E 05 0®5?904E 00 0®839l8E 02 0o53711E-05 

— 0® 20335E-03 0-26799E 00 0®21976E-01 0.11861E 00 

0.11I15E 01 0.46140E-01 0.58375E 00 0®75809E 00 

0-13486E 01 0.98508E 00 0.79679E-01 0«I0158E 00 

0-36130E 04 0.61566E 07 0®l7491E-04 0-60883E-06 

0®23024E 01 0 ® 10380E 06 0.61827E-03 0«64175E 02 

0-14235E 05 0®58557E 00 0®85481E 02 0®81709E-06 

— 0® 11 120E— 03 0®26798E 00 0. 19366E-01 0.11325E 00 

n 11600E 01 0® 14658E— 01 0®59837E 00 0.78306E 00 

® 14082E 01 0® 11238E 01 0®70536E-01 0.91882E-01 

0o36709E 04 0.59459E07 0- 16033E-04 0.59489E-06 

0.23803E 01 0.98931E 05 0®66338E-03 0.65629E 02 

0® 14738E 05 0.59855E 00 0.88597E 02 -0. 10283E-05 

0. 11303E— 03 0.26797E 00 0.14728E-01 O-1O282E0O 

0.12559E 01 — 0.52020E-01 0.62207E 00 0.83071E 00 

0® 15274E 01 0.I4736E 01 0.54089E-01 0.73863E-01 
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0.37874E 04 0.55116E 07 0.13263E-04 0.56531 £-06 

0.255Q8E 01 0.88858E-05 0.77322E-03 Q.687G7E 02 

0.156716 05 0.62424E 00 0.94780E 02 C.12220E-04 

— 0.58 288 E— 04 0.26795E 00 0. 726306-02 0.83295E-G! 

0.13506E 01 -0.12257E 00 0.63796E 00 C.87523E 00 

0* 16466E 01 0.I9588E 01 G.40285E-01 0.57947E-01 

0.39031E G4 0.50666E 07 0.10746E-04 0.53367E-06 

0.27417E 01 0.78592E 05 0.91725E-G3 0.720S9E 02 

0.165Q3E 05 0.64963E 00 0.10088E 03 0.68483E-05 

— 0.32718E— 03 0.26795EOO 0.14612E-02 0.65869E-01 

0.13976E 01 — 0.I5890E 00 0.642 87 E 00 0.S9627E 0© 

0.17062E 01 0.22685E 01 0.34391E-01 0.5G839E-G1. 

0.396G0E 04 0.48427E 07 0.95978E-C5 0.5I719E-06 

0.284536 01 D.73489E 05 0.10060E-02 0.73930E 02 

0.16881E 05 0.66222E 00 0.10388E 03 -Q.17S28E-04 

— Q . 9 1 44 0 E— 0 4 C.26794E 00 -0.94549E-03 0.58022E-0i 

0. 149 IDE 01 -0.23286E 00 0.64658E 00 0.93562E 00 

0.18254E 01 0.3O422E 01 0.24588E-01 0.38417E-01 

0 .407036 04 0.44G00E 07 0.75524E-05 0.48337E-06 

6.3Q68IE 01 0.63596E 05 0.12263E-02 Q.77988E 02 

0. 115616 05 0.68723E 00 0.10981E 03 -0.30811E-04 

0 .284136—03 0.26794E 00 -0.46642E-02 0.441986-01 

0.158376 01 -G.30755E 00 0.64225E 00 0.97059E GO 

0. 19446E. 01 0.39704E 01 0.I7299E-0I 0.2S481E-0I 

0 .417236 04 0.39794E 07 0.58750E-05 0.44956E-06 

G.33G70E 01 0.54526E 05 0.15148E-C2 0.82594E 02 

0.1S147E 05 0.71209E 00 0.11559E 03 0.36463E-04 

0. 2.1228c— 03 G.26793E GO -0.71896E-02 Q.33019E-01 

O.i 6762E 01 -G.38180E 00 0.63023E 00 0.99979E 00 

0.20638601 0.47943E 01 0.12253E-01 0.21054E-01 

C.42609E 04 0.36060E 07 0.45922E-05 0.41796E-06 

0. 354786.' 01 0.46815E 05 0.18719E-02 0.87634E 02 

C.18&49E 05 0.73686E 00 0.12I21E03 0.31856E-04 

-0.1.9119E-04 0.26792E 00 -0.88344E-02 0.24576E-07 


0. 17224E 01 -0 

0.21234E 01 
0.42983E 04 
0.36610E 01 
0.18873E 05 
-0.60526E-04 


-41841E 00 
0.49506E 01 
0.34461E 07 
0.43628E 05 
0.74924E 00 
0.26792E 00 


0.62153E 00 
0.10455E-01 
0.41001E-05 
0.20672E-02 
0.12395E 03 
-0.92038E-02 


0.1IO117E 01 
0« 18229E— 01 
0.40394E-06 
0.90187E 02 
0.1 1001 E-04 
0.21342E— 01 


0.18150E 01 — 0.48999E 00 0.59925E 00 0.10286E 01 

0.22426E 01 0.40800E 01 0.80854E-02 0.14170E-01 

0.43546E 04 0.32022E 07 0.34125E-05 0.38190E-06 

0.38477E 01 0.38912E 05 0.24369E-02 0.94824E 02 

0.19281E 05 0.77407E 00 J 0.12927E 03 0.4Q367E-05 

— 0.902 37E— 04 0.26793E 00 -0. 89484E-02 0.16665E-0I 


0.18&16E 01 — G.52472E 00 0.58594E 00 0.1Q330E 01 

0.23022E 01 0.29546E 01 0.74264E-02 0.12800E-01 

0.43723E 04 0.31253E 07 0.321I4E-05 0.37478E-06 

0.39105E 01 0.37465E 05 0-25798E-02 0.96649E 02 

0.1947 IE 05 0.78655E 00 0.13185E 03 0.79463E-05 

— 0.10557E— 03 0.26793E 00 -0.84302E-02 0.15075E-01 


0.19555E 01 -0.59162E 00 0.55575E 00 0. 10336E 01 

0.24214E 01 0.56508E 00 0.68950E-02 0.10809E-01 

0.43872E 04 0.30597E 07 0.30455E-05 0.36865E-06 

0.39657E 01 0.36245E 05 0.27321E-02 0.99026E 02 

0.19835E 05 0.81171E 00 0.13679E 03 0.38358E-05 

-0.14450E-03 0.26794E 00 -0»68359E-02 0.12745E-01 


0.20507E 01 -0.65476E 00 0.52175E 00 0.10276E 01 

0.25406E 01 0.0 0.68858E-02 0.98059E-02 

0.43875E 04 0.30586E 07 0.30426E-05 0.36854E-06 

0.39667E 01 0.36 223E 05 0.2772ZE-02 0.10042E 03 

0.20195E 05 0.83722E 00 0.1J4145E 03 -0.47177E-05 

— 0 . 62362E— 04 0.26794E 00 -0.47573E-02 0.11563E-01 


0.21472E 01 
0.26598E 01 
0.43867E 04 


-0.71409E 00 
0.10000E-02 
0.3062 IE 07 


0.48475E 00 
0.69134E— 02 
0.30513E— 05 


0.10212E 01 
0. 97204E— 02 
0. 36887E— 06 
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0.39637E 01 
0.2O555E 05 
0.42135E-04 


0.36288E 05 
0 o 86308E 00 
0.26794E 00 


0.28038E-Q2 
0. 14583E 03 
-0. 29891 E-02 


0.10174E 03 
-0.80277E-05 
0. 11462E— 03 


0.23438E 01 -0.82193E 00 0.40372E 00 O.I0096E Oi 

0. 28982E 01 0.I5596E-01 0.68900E-02 C.95088E-02 

0.43874E 04 G.30591E 07 0.30439E-05 0„36859E~‘06 

0- 39663E 01. 0.36233E 05 0.28849E-02 0.10453E 03 

0.21275E 05 G.91573E 00 0.15384E 03 -0„ 1G058E-0. 

-0.25259E-04 0.26794E 00 -0.36115E-03 0.11213E-DI 

0® 2682 16 01 -0.97376E 00 0.25415E GO C.93916E 00 

G.329S2E 01 0.2652GE-02 0.68941E-02 0..90523E-02 

G.43873E 04 G .305966 07 0.30452E-G5 0.36864E-06 

0-39658E 01 0.36243E 05 0.30122E-02 C.10917E 03 

0.22433E 05 0.10064E 01 0.16537E 03 0.23836E-04 

-0.18715E-03 0.26796E 00 0.17401E-02 0.10674E-01 


0.3G289E 01 — 0.I0951E 01 0.96266E-01 0.96890E GO 

0.369S2E- ‘01 0*0 0.6892IE-02 G.86548E-02 

0.43873E 04 0.30594E 07 G.3G446E-05 0.3686IE-06 

0.39660E 01 0. 36238E 05 0.31418E-02 0.11385E 03 

0*2369 IE C?5 0.10993E 01 0.17498E03 -0.59520E-04 

0® 13732 £-03 0-26795E0G 0.11290E-02 0.10206E-0I 

0.32048E 01 -Q.U463E 01 0.I6196E-01 0.95875E 00 

0.3E982E 01 0*0 0.68923E-G2 0.84793E-02 

0.438T3E 04 G.30594E 07 C.30447E-05 0.36861E-06 

0.39660E Oi 0.36238E 05 0.32068E-02 0.1162IE 03 

G» 24295E 05 0.I1465E 01 0.17919E 03 -G.15637E-04 

0®31 380E— 03 0.26795E 00 0.47901E-03 0.99987E-02 

0.32269E 01 -G„il524E 01 0.61737E-Q2 0. 95748E 00 

0* 39232c 01 0.0 0.68923E-02 0.84578E-02 

G.43873E 04 0.30594E 07 0.30447E-05 0.3686i£H)6 

0.39660E 01 0.36238E 05 0.32150E-02 0.11651E 03 

0.24371E 05 G.11524E 01 0.17969E 03 -Q.45516E-C5 

0.3I7I8E-03 0.26795E 00 G.41198E-G3 0.99733E-02 
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32380E 01 — 0.11553E 01 0.11631E-02 0.95684E 00 

-39357E 01 0.12786E— 04 0.68923E-02 0.84476E-02 

0.43873E 04 0.30594E 07 0.30447E-05 0.36861 E-06 

0.39660E 01 0.36238E 05 0.32I91E-02 0.11665E 03 

0.24408E 05 0.11553E 01 0.17994E 03 -0. 11849E-05 

0.31737E-03 0.26795E 00 0.37953E-03 0.996I4E-02 

0.32408E 01 -0.11561E 01 0.L6161E-04 0.95669E 00 

0.39388E 01 0.47180E-04 0.68923E-02 0. 84458E.-02 

0.43873E 04 0.30594E 07 0. 30447E-05 0.36861E-06 

0.39660E 01 0.36238E, 05 0.32201E-02 0.11669E 03 

0.24418E 05 0.11561E 01 0.18000E 03 0. 18591E-06 

0.31727E-03 0.26795E 00 0.37152E-03 0.99592E-02 

0.32421E 01 -0.11564E 01 -0.14044E-03 0.95661E 00 ; 

0.39404E 01 0.43983E-04 0.68923E-02 0.84444E-02 

0.43873E 04 Q.30594E07 0.30447E-05 0.36861E-06 

0.39660E 01 0.36238E 05 0-32206E-02 0.11671E 03 

0.24423E 05 0.11564E 01 0.18001E 03 0.37077E-06 

0.31720E— 03 0.26795E 00 0.36762E-03 0.99575E-02 

0. 32449 E 01 -0.11572E 01 -0.14053E-03 0.95645E 00 

0.39435E 01 0.42898E-04 0.68923E-02 0.84417E-02 

0.43873E 04 0.30594E 07 i 0.30447E-05 0.36861 E-06 

0.39660E 01 0.36238E 05 0.32216E-02 0.11675E 03 

0.24432E 05 0.11572E 01 0.18001E 03 0.36192E-06 

0.31707E-03 0.26795E 00 0.35966E-03 0.99544E-02 

0.32504E 01 -0.11587E 01 -0.14071E-03 0.95613E 00 : 

0.39497E 01 0.40743E-04 0.68923E-02 0.84363E-02 

0.43873E 04 0.30594E 07 ^ 0.30447E-05 0.36861 E-06 

0.39660E 01 0.36238E05 0.32237E-02 0.11682E 03 

0.24451E 05 0.11587E 01 0.18001E 03 0.28289E— 05 

0.31676E— 03 0.26795EOO 0.34392E-03 0.99480E-02 

0.326I5E 01 — 0-11616E 01 -0.14107E-03 0.95550E 00 

0.39622E 01 0.36502E-04 0.68923E-02 0.84255E-02 

0.43873E 04 0.30594E 07 0.30447E-05 0.36861E-06 

0.39660E 01 ; 0.36238E 05 0.32278E— 02 0.11697E 03 
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Q.24489E 05 
0-31604E— 03 


0.11616E 01 
G.26795E 00 


G.18001E 03 
0.313G2E-03 


0. 27622E-05 
Q.99353E-02 


0-32836E 01 -0.1I676E01 -0.14179E-03 0.95423EOO 

0.39872E OS 0.28294E-04 0.68923E-02 0.8404QE-Q2 

0-43873E 04 0.30594E 0? 0.30447E-05 G.36861E-0& 

G.3966GE 01 0.36238E 05 0.32359E-02 Q.11727E 03 

G.24564E 05 0.11676E 01 C.1800IE 03 0.94361E-0? 

0.314I6E— ©3 0.26795E 00 0.25341E-G3 0.99099E-G2 

0.33279E 01 -0.11794E 01 -C.14323E-03 0.95169E 00 

0.40372E 01 0* 12957E—04 0.68923E-02 G.83613E-Q2 

0.43873E 04 0.30594E 0? 0.30447E-05 ©. 36861 E-Q6 

0.39660E 01 0.36238E 05 0.32523E-Q2 0-11786E 03 

0.24715E 05 0.11794E 01 0.18001E 03 0.19234E-05 

0.3G873E— 03 0.26795E 00 0.14289E-03 0.98595E-02 

S.34163E 01 — 0® 12031E 01 -0.14611E-03 0.94663E 00 

0.41372E 01 0.0 0.68923E-02 0.82772E-02 

0.43873E 04 0.30594E 07 C.30447E-05 0.36861E-06 

0.39660E 01 0.36238E 05 0.32851E-C2 0.11905E 03 

0.25017E 05 0.12031E 01 0.18001E 03 -0.24624E-G6 

0.29186E-03 0.26795E 00 -0.44488E-04 0.97603E-02 


0.35933c 01 -0.125Q5E 01 -0.15187E-03 0*9365 IE 00 

0.43372E 01 0.0 0.68923E-02 0.81143E-02 

0.43873E 04 0.3G594EG7 0.30447E-05 0. 36861 E-06 

0.39660E 01 0.36238E 05 0.3351IE-G2 0.12144E 03 

0.256 2 IE 05 0.12505E 01 0.18001E 03 -0.52000E-06 

0.238896-03 0.26795E00 -0.29365E-03 0.95683E-02 

0.3947 BE 01 -0.13454E 01 -0.16339E-03 0.91633E 00 

0-47372E 01 0.0 0.68924E-02 0.78039E-02 

0.43873E 04 0.30594E 07 0.30447E-05 0.36861E-06 

0.39660E 01 0.36238E 05 0.34844E-02 0.12627E 03 

0.26829E 05 0.13454E 01 0.18001E 03 -0.52600E-06 

0.893426-04 0.26795E 00 -0.35161E-03 0.92023E-02 

0.43012E 01 -0.I4402E 01 -0.17491E-03 0.89621E 00 
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0.51372E 01 0.0 0.68924E-02 0.75119E-02 

0.43873E 04 0.305 94E 07 0.30447E-05 0.36861 E-06 

0.39660E 01 0.36239E 05 0.36198E-02 0.13118E 03 

0.28037E 05 0.14402E 01 0.18001E 03 -0.27437E-05 

-0.68888E-04 0.26795E 00 -0.11996E-03 0.88579E-02 

0.46552E 01 -0.15351E 01 -0.18642E-03 0.87614E 00 

0.55372E 01 0.7I050E-05 0.68924E-02 0.72363E— 02 

0.43873E 04 0.30594E 07 0.30447E-05 0.36861 E-06 

0.39660E 01 0.36239E 05 0.37579E-02 0.136I8E 03 

0.29245E 05 0.15351E 01 O.IBOOIE 03 -0.15874E-Q5 

-0.20417E-03 0.26795E 00 0.22132E-04 0.85329E-02 

0.50092E 01 — 0.16299E 01 -0.19794E-03 0.85612E 00 

0.59372E 01 0.19909E-04 0.68924E-02 0.69749E-02 

0.43873E 04 D.3G594E 07 G.30447E-05 0.36861 E-06 

0.39660E 01 G.36239E 05 0.38989E-02 0.14129E 03 

0.30453E 05 0.16299E 01 0.1800IE 03 0. 13830E-05 

— 0.29784E— 03 0.26795E 00 0.84565E-04 0.82247E-02 

0.53631E 01 -0.17248E 01 -0.20946E-03 0.83614E 00 

0.63372E 01 3.19395E-04 0.68924E— 02 0.67257E-02 

0.43873E 04 0.30594E07 0. 30447E-05 0.36861 E-06 

0.39660E 01 0.36239E- 05 0.40433E-02 0.14652E 03 

0.31661E 05 0.17248E 01 0.18001E 03 0.1 6061 E-05 

-0.33534E-03 0.26795E 00 0.80416E-04 0.79309E-02 

n 57171E 01 -0.18196E 01 -0.22098E-03 0.81619E 00 

.67372E 01 0.85929E-05 0.68924E-02 0.64875E-02 

0.43873E 04 0.30594E 07 0.30447E-05 0.36861E-06 

0.39660E 01 -- 0.36239E 05 0.41916E-02 0oi5190E 03 

0.32869E 05 0.18196E 01 0.18001E 03 0. 14535E-05 

— 0.30688E-03 0.26795E 00 0.33363E-04 0.76500E-02 

0.6071 IE 01 -0.19145E 01 -0.23250E-03 0.79626E 00 : 

0.71372E 01 0.0 0.68924E-02 0.62593E-02 

0.43873E 04 0.30594E 07 0.30447E-05 0.36861 E-06 

0.39660E 01. 0.36238E 05 0.43442E-02 0.15743E 03 

0.34077E 05 0.19145E 01 0.1 800 IE 03 0.10148E-05 
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-G.22756E-03 C.26795E 00 -0.560416-06 0.73809E-02 

O.64250E ©1 -G.20G93E 01 -G.24401E-03 CUT7637E00 

0.753725 01 0.0 0.68924E-02 0.60405E-G2 

0.438736 D4 0-3G594E 07 0.30447E-05 0.368616-06 

0.39660E 01 0.36239E 05 0.45016E-02 0.163136 03 

G.352S5E 05 0.20093E 01 Q.18001E 03 -0.18226E-G6 

-0. 125926-03 0.26794E 00 -0.183775-04 0.712296-02 

0.67790E 01 — 0.2 1041 E 01 -0.25553E-03 0.75650E 00 

0.79372E 01 0.0 0.68924E-02 G.58298E-02 

0.438736 04 0.30594E 0? 0.30447E-05 G.36861E-0© 

O.39660E €1 0.36239E 05 Q.46643E-0? . 0.16903E 03 

0.364936 05 0.21041E 01 0.18001E 03 0.13815E-06 ' 

— 0.25706E— 04 0.26795E 00 -0.22336E-04 0.68745E-02 

0.71329E 01 -0.219906 01 -0.26705E-03 0.736&5E 00 

0.S3372E 01 O.O 0.68924E-02 0.56264E-02 

0.43873E 04 0-30594E 07 O.3O447E-05 0.36861E-Q6 

0. 396605 Cl 0.36239E 05 0.48329E-C2 0.I7514E 03 

0.37701E 05 0.2I990E 01 0.18001E 03 -G.27423E-0? 

0.-S3610E-C4 0.26795E 00 -0.15003E-04 0.66346E-02 

0.74869E 01 -0.22938E01 -0.27857E-03 O.T1682EOO 

0.87372E 01 0.0 0.68924E-02 0.54296E-02 

0.438735 04 0.30594E 07 0.30447E-05 0.368615-06 

0.39660c Cl 0.36239E 05 0.50080E-02 G.18148E 03 

G.3S90SE 05 0.229385 01 0.18001E 03 -0.30999E-06 

0.104195-03 0.26795E 00 -0.71848E-05 Q. 640266-02 

0.78409E 01 — 0.23887E 01 -0.29009E-03 G.69700E00 

0.91372E 01 0.0 0.68924E-02 0.52389E-02 

0.43873E 04 0.30594E 07 0.30447E-05 G.36861E-06 

0.39660E 01 0.36239E 05 0.51904E-02 0.18809E 03 

0.40H6E 05 G.23887E 01 0.18001E 03 -0.97911E-07 

0.130635-03 0.26795E 00 -0.11603E-05 0.61776E-02 

0.81946E 01 -0.24835E 01 -0.30160E-03 0.67720E 00 

0.95372E 01 0.55306E— 06 0.68924E-02 0.50536E-02 
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0.438736 04 Q.30594E 07 0o30447E-05 0.368616-06 

0.39660601 0.36239E 05 0.53807E-02 0.19499E03 

0.413246 05 0.248356 01 0.180016 03 0.15816E-06 

0.13760E-03 0.267956 00 0.29880E-05 0.595916-02 

0.854886 01 -0.257846 01 -0.313126-03 0.65741E 00 

0.99372E 01 0.17939E-05 0.689246-02 0.48732E-02 

0.43873E 04 0.30594E 07 0.30447E-05 0.36 861 E-06 

0-39660601 0.36239E 05 0.55799E-02 0.2022 IE 03 

0. 42532 E 05 0.25784E 01 0.18001E 03 0. 11922E-05 

0.12910E-03 0.26795E 00 0-56868E-05 0.57464E-02 


0.89028E 01 — 0.26732E 01 -0. 32464E-03 0.63764E 00 

0.1Q337E 02 ‘ 0.26430E-05 0.68924E-02 0.46973E-D2 

0.43873E 04 0.30594E 07 0. 30447E-05 0.36861E-06 

0.39660t 01 0.36239E 05 0.57889E-02 0.20978E03 

0.43740E 05 0.26732E 01 0.1800IE 03 -0.28392E-05 

0.10855E-03 0.26795E 00 0.69790E-05 0.5539OE-02 

0.92567E 01 -0.27681E0I -0.33616E-03 0.61787E 00 

0.10737E 02 9.30621E-05 0.68924E-02 0.45255E-02 

0.43873E 04 0.30594E 07 0. 30447E-05 0.36861E-06 

0-39660E 01 0.36239E 05 0.60088E-02 0.21775E 03 

0.44948E 05 0.27681E 01 0.18001E03 0.28173E-05 

0.78905E-04 0.26795E 00 0.70305E-05 0.53364E-02 

0.96107E 01 — 0.28629E 01 -0.34768E-03 0.59812E 00 

0.UI37E 02 0. 30038E— 05 0.68924E-02 0.43573E-02 

0.43873E 04 0.30594E 07 0.30447E-05 Q.3686IE-06 

0.39660E 01 0.36239E 05 0.62406E-02 0.22615E 03 

0.46156E 05 0.28629E 01 0.18001E 03 -0.268146-06 

0.42735E-04 0.26795E 00 0.595446-05 0.51381E-02 

° 996466 01 -0.295776 01 -0.35920E-03 0.57837E 00 

.115376 02 0.24096E— 05 0.68924E-02 0.41926E-02 

0.43873E 04 0.30594E 07 0.30447E-05 0.368616-06 

0.39660E 01 0.36239E 05 0.648586-02 0.23504E 03 

0.47364E 05 0.29577601 0.18001E 03 -0. 18351E-05 

0.22725E-05 0.26795E 00 0.38591E-05 0^49439E-02 
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0.10319E 02 -G.30526E 01 -0.37071E-03 0. 

0- 11937E 02 0. 150526-05 0.68924E-02 € 

G.43873E 04 0.3C594E 07 0.30447E-05 

0.39660E 01 0.36239E 05 0.67459E-02 

0.48572E 05 0.30526E 01 0-18001E 03 

-0.40459E-04 0.26795E 00 0.14742E-05 

0.1O673E 02 — G.31474E 01 -0.38223E-03 0< 

0.12337E 02 0.69834E-06 0.68924E-02 C 

0-43873E 04 0.30594E 07 0.30447E-05 

0.39660E 01 0.36239E 05 0.70227E-02 

0.49780E 05 0.31474E 01 0.18001E 03 

-0.83332E-04 0-26795E 00 -G.34088E-06 

O® 1 1027E 02 — Q.32423E 01 -0.39375E-03 0. 

0.12737E 02 G.63295E-08 ©.68924E-02 ( 

0.43873E 04 0.30594E 07 0.30447E-05 

0.39660E 01 Q.36239E 05 0.73180E-02 

0-50988E 05 0.32423E 01 0.18001E 03 

-0.12435E-03 0.26795E 00 -0.17005E-05 

0.I138GE 02 -0®33371E 01 -0.40527E-03 0. 

0® 13137E 02 OoO 0.68924E-02 « 

0.43873E 04 Q®30594E 07 0.30447E-05 

0.39660E 01 0.36239E 05 0.76344E-02 

0.52196E 05 0.33371E 01 0-18001E 03 

—0.151752—03 0.26795E 00 -0.26166E-05 

0.11734E 02 -Q.34320E 01 -0.41679E-03 0. 

0.13537E 02 0.0 0.68924E-02 I 

0.43873E 04 0.30594E 07 0. 30447E-05 

0.39660E 01 0-36239E 05 0.79744E-02 

0.53404E 05 0.34320E 01 0.18001E 03 

— 0® I9396E— 03 0.26795E 00 -0.30567E-05 

0.1 2088 E 02 -0.35268E 01 -0.42830E-03 0. 

0.13937E 02 0.0 0.68924E-02 i 

0.43873E 04 0.30594E 07 0.30447E-05 


55863E 00 
<®40309E— 02 
0® 36881 E-06 
0.24446E 03 
0.467 17 E-06 
0.47532E-02 

53890E 00 
U38720E— 02 
0.36861 E-06 
0.25449E 03 
0. 12640E— 05 
0.45659E-02 

• 51918E 00 
3.37157E— 02 
0.368616-^-06 
0.26519E 03 
0.1 9998 E-05 
0.43815E— 02 

.49946E 00 
3.35617E— 02 
0.36861 E-06 
0.27666E 03 
0. 22260E— 05 
0.42000E-02 

.479746 00 
Q.34099E— 02 
0.36861E— 06 
0-28898E 03 
-0.26574E-05 
0.40209E-02 

. 46004E 00 
0.32599E-02 
0.36861E— 06 
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0.39660E 01 0.36239E 05 0.83412E-02 0.30227E 03 

0o546l2E 05 0.35268E 01 0.18001E 03 -0.12199E-05 

-0.21958E-03 0.26795E 00 -0.31604E-05 0.38441E-02 

0.12442E 02 -0.36217E 01 -0.43982E-03 0.44033E 00 

0. 14337E 02 0.0 0.68924E-02 0.31117E-02 

0.43873E 04 0.30594E 07 0.30447E-05 0.36861E-06 

0.39660E 01 0.36239E 05 0.87385E-02 0.31667E 03 

0.5582QE 05 0.36217E 01 0.18001E 03 0.24390E-OS 

7^0 • 2 3735E-0 3 Q.26795E 00 -0.28554E-05 0.36693E-02 

“ 12796E 02, -0.37165E, 01 -0.45134E-03 0.42063E 00 

. 14737E 02 0.0 0-68924E-02 0.29650E-02 

0.43873E 04 0.30594E 07 0.30447E-05 0.3686TE-06 

0.39660E 01 0.36239E 05 0.91708E-02 0.33234E03 

0.57028E 05 0.37165E 01. 0.18001^ 03 0.75820E-06 

-0.24615E-03 0.26795E 00 -0.22694E-05 0.34963E-02 

13150E 02 -0.38113E Ol -0.46286E-03 0.40094E 00 

. 13 137 E 02 0.0 ■ 0. 68924E— 02 0.28198E-02 

0.43873E 04 0.30594E 07 0.30447E-05 G.3&86IE-06 

0.39660E 01 0.36239E 05 0.96432E-02 0.34946E03 

0.58236E 05 0.38113E 01 0.18001E 03 -0.24924E-05 

-0.24500E-03 /0.26795E 00 -0.14388E-05 0.33251E-02 

13504E 02 -0.39062E 01, -0.47438E— 03 0.38124E 00 

. 15537E 02 0.150G3E— 06 0.68924E-02 0.26758E-02 

0.43873E 04 0.30594E 07 0.30447E-05 0.36861 E-06 

0.39660E 01 0.36239E 05 O.L0162E-O1 0.36826E03 

0.59444E 05 0.39062E 01 0.18001E 03 -0.14432E-05 

— 0 .2 3445E— 0 3 0.26795E 00 -0. 66377E-06 0.31553E-02 

Q.13858E 02 -0.40010E 01 -0.48589E-03 0.36155E 00 

0.15937E 02 0.53792E-06 0.68924E-02 0.25329E-02 

0.43873E 04 0.30594E 07 0.30447E-05 0.36861 E-06 

0.39660E 01 0.36239E 05 0.10735E-01 0.38903E03 

0.60652E 05 0.40010E 01 0.18001E 03 -0.21730E-05 

-0.21669E-03 0.26795E 00 -0.59589E-07 0.29868EM>2 
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0.14212E 02 -0.40959E 01 -0.49741E-Q3 

0.16337E 02 0.75467E-06 0.68924E-02 

0.43873E 04 0.30594E 07 0.30447E-05 

0.39660E 01 0.36239E 05 0.11372E-01 


0.34187E 00 
0.23911E— 02 
0. 36861 E~06 
0.4121 IE 03 


0.61859E 05 
— G. 19383E— 03 


G.40959E 01 
0.26795E 00 


0.18001E 03 
0. 42930E— 06 


-0. 38948 E-05 
0.28195E— 02 


0.14566E 02 -0.41907E 01 -0.50893E-03 0.32218E 00 

0.16737E 02 0.74550E-06 0.68924E-02 0«22501E-02 

0.438736 04 0.30594E 07 0.30447E-05 0.36861E-06 

0.39660E 01 0.36239E 05 0.12085E-01 0.43794E 03 

0.63067E 05 0.41907E 01 0.18001E 03 -0.52045E-05 

— 0. 16775E— 03 0.26795E 00 Q.81071E-06 0.26533E-02 

0.1492OE 02 — 0.42856E 01 -0.52045E-03 0.30250E 00 

0.I713TE 02 0.43370E— 06 0.68924E-02 0o2I099E-02 

0.438 T3E 04 Q.30594E 07 0.30447E-05 0-36861 E-06 

0.39660E 01 0.36239E 05 0.12888E-01: 0.46703E 03 

0.64275E 05 0.42856E 01 0.18001E 03 -0.68415E-05 

— 0. 14020E— 03 0.26795E 00 0.11127E-05 0.24880E-02 


0.15274E 02 -0.43804E 01 -0.53196E-03 0.28282E 00 

0.17537E 02 0.0 0-68924E-0? 0.19704E-02 

0. 43873 E 04 0.30594E 07 0.30447E-05 0.36861E-06 

0.39660E 01; 0.36239E 05 0.13800E-01 0.500IGE 03 

0.65483E 05 0.43804E 01 0.18001E 03 -Q.94807E-05 

— 0. 112T5E— 03 0.26795E 00 0.12383E-05 0.23235E-02 

0.15628E 02 -0.44752E01 -0.54348E-03 0.26314E00 

0.17937E 02 0.0 0-68924E-02 0.18315E-02 

0.43873E 04 0.30594E 07 0.30447E-05 0.36861E-06 

0.39660E 01 0.36239E 05 0.14847E-01 0.53802E 03 

0. 66691 E 05 0.44752E 01 0.18001E 03 0.68809E-05 

— 0. 86860E— 04 0.26795E 00 0.13652E-05 0.21597E-02 

0.15982E 02 -0.45701E 01 -0.5550DE-03 0. 24346E 00 

0.18337E 02 0.0 0.68924E-02 0.I6931E-02 

0.43873E 04 0.30594E 07 0.30447E-05 0.36861 E-06 

0.39660E 01 0.36239E 05 0.16060E-01 0.58200E 03 
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0©67S99E 05 
-0.63853E-04 


0.45701E 01 
0.26795E 00 


0.18001E 03 
0.13434E-05 


— 0.10090E— 04 
Q.19965E— 02 


*** SUBSCRIPT OUT OF RANGE *** 

CALCULATIONS FOR STREAMLINE NO. 2 DISCONTINUED 

THE CALCULATIONS FOR 2 STREAMLINES HAS BEEN COMPLETED 
PROGRAM TERMINATING 
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